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ABSTRACT 


An  aerodynamic  stability  computer  program  (SFFSTAB)  has  been 
developed  for  calculating  the  spin  rate  required  for  stabilization  of  a 
self-forging  fragment.  An  aerodynamic  stability  criterion  which 
combined  gyroscopic  and  dynamic  stability  was  used  together  with  a 
technique  for  calculating  aerodynamic  coefficients.  SFFSTAB  is  a 
useful  tool  for  conducting  aerodynamic  stability  parameter  studies  for 
different  fragment  shapes.  Complete  documentation  of  the  computer 
program  including  sample  problem,  flowchart  and  FORTRAN  listing  is 
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P 

R 

R1 
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projectile  cross-sectional  area 


dynamic  stability  factor 


gyroscopic  stability  factor 

velocity 
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1.0  INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  criteria  for 
stability  and  to  document  a  program  which  is  capable  of  calculating 
aerodynamic  coefficients  and  spin  rates  required  for  stabilization  of 
self-forging  fragments. 

The  Demolitions  Group  at  the  Defence  Research  Establishment 
Suf field  (DRES)  is  responsible  for  military  related  research  in  Canada 
on  demolitions  and  demolition  devices.  DRES  is  developing  a 
Self-Forging  Fragment  (SFF)  technology  base  suitable  for  the  study  of 
the  application  of  SFF  devices  to  demolition  problems. 

One  of  the  projects  currently  under  stucfy  is  the  application  of 
spin  stabilization  to  self-forging  fragments.  The  project  requires 
both  numerical  and  experimental  capabilities  to  study  the  SFF  formation 
and  flight  stability. 
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The  ORES  Flash  X-ray  Facility  has  the  capability  of  obtaining 
formation,  shape,  stability  characteristics  and  trajectory  data  for  the 
fragments.  The  fragments  can  be  monitored  over  a  60  m  flight  path  by 
using  a  combination  of  X-ray  radiographs  and  yaw  screens.  To  fjlfill 
the  numerical  requirements,  and  the  requirement  for  a  code  to  determine 
what  parameters  to  vary,  a  computer  program  was  developed  to  calculate 
the  aerodynamic  coefficients  for  a  self-forging  fragment  and  to  predict 
its  stability. 

Section  2  describes  self-forging  devices  and  the  fragment 
characteristics,  including  selected  results  from  experimental  tests. 

In  Section  3  aerodynamic  stability  criteria  are  presented.  Section  4 
describes  the  computer  program  Self-Forging  Fragment  STABility 
(SFFSTAB).  The  results  for  sample  problems  are  also  described  in  this 
section. 


The  DATCOM  [3]  method  for  obtaining  aerodynamic  coefficients  is 
documented  in  Appendix  A,  and  properties  of  spherical  segments  and 
truncated  cones  are  given  in  Appendix  B.  A  flowchart  and  a  FORTRAN 
listing  of  the  computer  program  with  input  and  output  for  the  sample 
problem  presented  in  Section  4  are  contained  in  Appendix  C. 


2.0  SELF-FORGING  FRAGMENT  DEVICE 

A  typical  self-forging  fragment  (SFF)  device,  as  illustrated  in 
Figure  1,  consists  of  a  body  which  is  filled  with  high  explosive  (HE) 
and  a  metal  liner.  The  charge  is  initiated  by  a  detonator  and  an 
explosive  booster.  Upon  initiation,  the  resulting  detonation  wave 
travels  through  the  HE  subjecting  the  metal  liner  to  high  forces 
causing  a  single  fragment  to  be  formed,  and  accelerating  this  fragment 
to  a  high  velocity.  Typical  SFF's  have  initial  velocities  of 
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i-i  kin/sec.  The  >e  fragments  ^an  travel  to  distances  well  in  excess  of 
lclOo  calibers,  beyond  which  accuracy  becomes  a  problem. 

In  many  sel f-f arging  fragment  designs,  the  objective  is  to 
produce  a  long  rod-type  fragment  as  opposed  to  a  ball -type  fragment 
since  the  former  has,  in  general,  greater  target  penetration  capabil¬ 
ity.  However,  a  long  fragment  is  inherently  unstable  in  flight  as 
illustrated  in  Figure  2.  This  figure  shows  flash  radiographs  of  the 
fragment  as  it  deviates  from  tne  intended  flight  path.  This  instabil¬ 
ity  results  in  a  loss  of  accuracy  and  a  lower  target  penetration  cap¬ 
ability.  Current  SFF  designs  use  a  flared  tail  fragment.  This  feature 
increases  the  stability  of  the  fragment  and  improves  its  accuracy. 
However,  the  flared  tail  increases  the  aerodynamic  drag  on  the  frag¬ 
ment,  resulting  in  a  shorter  flight  path,  and  it  also  reduces  the  pene¬ 
trating  capability  of  the  fragment  since  the  tail  does  not  contribute 
to  the  penetration  process.  If  a  fragment  could  be  spin  stabilized, 
the  requirement  for  a  flared  tail  would  be  eliminated  (see  Figure  3). 
This  would  result  in  increased  performance  of  the  self-forging  fragment 
device.  The  main  aim  of  this  report  is  to  present  a  computer  program 
for  calculating  the  requirements  for  aerodynamic  stability  and  to 
illustrate  its  application  to  self -forging  fragments  with  flared  rod 
and  blunt  cone  shapes. 


3.0  AERODYNAMICS  OF  SELF-FORGING  FRAGMENTS 


Determination  of  the  aerodynamic  stability  of  a  projectile 
requires  a  knowledge  of  the  mass  distribution  within  the  projectile 
( i . e . ,  location  of  the  center  of  gravity,  moments  of  inertia)  and  the 
aerodynamic  parameters  (i.e.,  lift,  drag,  etc.).  For  most  projectiles 
the  mass  distribution  can  be  calculated  or  measured  from  the  projectile 
design  and  the  aerodynamic  parameters  can  be  determined  by  wind  tunnel 
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experiments.  However,  for  self-forging  fragments  the  actual  shape  of 
the  fragment  that  will  be  formed  from  a  liner  is  unknown,  so  that  these 
properties  cannot  be  determined  a  priori.  Flash  X-radiography  can  be 
used  to  determine  the  outer  profile  of  the  fragment  but  in  general, 
unless  very  high-energy  X-rays  are  used,  the  radiographs  do  not  reveal 
any  information  about  the  interior  of  the  fragment.  Also,  the  high 
velocities  of  self-forging  fragments  requires  that  very  specialized 
hypersonic  tunnels  be  used  to  measure  the  aerodynamic  parameters. 
Numerical  techniques  are  therefore  used  to  calculate  the  self-forging 
fragment  stability.  This  report  describes  the  SFFSTAB  code  developed 
for  this  purpose.  After  a  review  of  the  aerodynamic  forces  and  moments 
acting  on  a  projectile,  formulas  for  gyroscopic  and  dynamic  stability 
used  in  a  stability  criterion  derived  by  Murphy  [2]  are  presented. 
Aerodynamic  coefficients  required  for  the  stability  formulas  are 
provided  by  the  DATCOM  method  [3].  Results  for  typical  fragment  shapes 
are  also  presented. 

The  aerodynamic  forces  acting  on  a  projectile  are  illustrated 
in  Figures  4  and  5.  The  resultant  (R)  of  the  forces  acting  on  a  Dro- 
jectile  can  be  decomposed  into  the  drag  force  (D)  and  lift  force  (L)  in 
the  wind  coordinate  system,  or  the  normal  force  (N)  and  axial  force  (A) 
in  the  body  coordinate  system.  The  resultant  force  (R)  acts  on  the 
center  of  pressure  (CP)  which  is  the  point  at  whicn  the  net  moment  is 
zero.  An  additional  force  which  acts  on  a  projectile  is  the  magnus 
force  (see  Figure  5).  The  combination  of  the  air  flowing  over  the 
projectile  and  the  spin  of  the  projectile  results  in  an  asymmetric 
boundary  layer  thickness  distribution.  This  effectively  changes  the 
aerodynamic  shape  of  the  body  which  results  in  a  force  normal  to  the 
angle  of  attack  plane. 

The  location  of  the  center  of  pressure  (CP)  depends  on  the 
velocity  of  the  projectile.  Therefore,  it  is  more  convenient  to  use 
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the  center  of  gravity  (CG)  or  center  of  mass  which  depends  only  on  the 
mass  distribution  in  the  projectile.  The  forces  acting  at  the  center 
of  pressure  are  replaced  by  an  equivalent  force  and  moment  combination 
at  tne  center  of  gravity.  The  aerodynamic  forces  and  moments  are 
usually  non-dimensional ized.  The  corresponding  force  and  moment  coef¬ 
ficients,  and  CmomenT.  take  the  general  forms: 


CF0RCE  =  F0RCE  7  !spV2S‘ 

(3.1) 

MOMENT  ~  MOMENT  /  *spV2Sl , 

(3.2) 

where  p  and  V  are  the  air  density  and  velocity  respecti vely,  S  is  the 
projectile  cross-sectional  area,  and  1  is  the  projectile  length. 

The  dominant  aerodynamic  effect  is  the  pitching-moment  or 
static-moment  coefficient  (Cm).  This  coefficient  represents  the  over¬ 
turning  or  pitching  moment  which  is  a  function  of  the  distance  between 
the  centers  of  mass  (CG)  and  pressure  (CP).  If  the  derivative  of  the 

pitching  moment  (C^  )  with  respect  to  the  angle  of  attack  (a)  is  nega- 

a 

tive,  the  CP  is  located  behind  the  CG  and  the  projectile  is  said  to  be 
statically  stable  (i.e.,  stable  as  long  as  it  is  not  perturbed).  Con¬ 
versely,  if  Cm  is  greater  than  zero,  then  the  CP  is  located  ahead  of 

a 

the  CG  and  the  projectile  is  statically  unstable.  Physically,  the 
center  of  gravity  is  the  pivot  point  of  the  projectile.  The  resultant 
force  (R)  acts  at  the  center  of  pressure  in  a  direction  towards  the 
rear  of  the  projectile.  Thus,  if  the  force  acts  ahead  of  the  pivot 
point  the  projectile  is  statically  unstable,  and  conversely  if  it  acts 
behind  the  pivot  point  the  projectile  is  statically  stable. 

Stabilization  of  projectiles  is  accomplished  by  two  techniques, 
fin  and  spin  stabi 1 ization.  Fin  stabil ization  is  accomplished  by 
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adding  fins  (or  a  flared  tail  in  the  case  of  a  self-forging  fragment) 
to  the  projectile  in  order  to  shift  the  center  of  pressure  behind  the 
center  of  gravity.  In  the  second  technique,  the  projectile  is  spun  in 
order  to  maintain  the  center  of  pressure  as  close  as  possible  to  the 
trajectory  which  is  the  path  of  the  center  of  mass.  For  a  spinning 
projectile  a  gyroscopic  stability  factor  (S^)  is  defined  by  Murphy  [2] 
as: 


S 


9 


2 


I  i 


x 


P  W  1  1 


(3.3) 


where  P  and  V  are  the  projectile  spin  rate  and  velocity,  is  the  air 
density,  d  is  the  projectile  diameter,  Cm  is  the  derivative  of  the 

Cl 

pitching-moment  coefficient,  and  I  and  I  are  the  axial  and  transverse 

x  y 

moments  of  inertia,  respectively.  The  dominant  aerodynamic  coefficient 
affecting  the  motion  of  a  spinning  projectile  is  the  pitching-moment 
coefficient.  In  order  for  a  statically  unstable  projectile  to  have  a 
periodic  motion  (i.e.,  continue  to  spin)  and  not  tumble,  the  gyroscopic 
stability  factor  (S^)  must  be  greater  than  one.  Physically  the  motion 
in  this  case  is  the  same  as  that  of  a  spinning  top.  Gyroscopic  stabil¬ 
ity,  like  static  stability  (i.e.,  center  of  pressure  behind  the  center 
of  gravity)  is  required  for  stable  oscillatory  motion,  but  this  alone 
does  not  ensure  that  small  perturbati ons  acting  on  the  projectile  will 
not  grow  and  result  in  unstable  motion.  In  order  to  characterize  the 
response  to  small  perturbations,  a  dynamic  stability  factor  (S^)  was 
derived  by  Murphy  [2]  from  the  equations  of  motion  of  a  spinning  pro¬ 
jectile.  This  stability  factor  is: 


I 
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Sd  = 


2  <C|(  -  CD)  +  2  Kx-2  Cm 

a  pu 

CM  -  2Cn  -  K  ‘MC  +  C  ] 
N  D  y  mq  m. 

u.  a 


(3.4) 


where  K  and  K  are  the  dimensionless  axial  and  transverse  radii  of 
x  y 

gyration,  respectively,  the  sum  (C  +  Cm  )  is  the  pitch-damping  coef- 

u 

ficient,  C..  and  C  are  the  derivatives  of  the  normal-force  and  mag- 
N  in  3 

a  p  - 

nus -moment  coefficients,  respectively,  and  CQ  is  the  drag  coefficient. 


The  relation  between  the  spin  rate  of  a  projectile  and  various 
types  of  projectile  stability  was  characterized  by  Murphy  [2].  His 
results  are  summarized  on  the  stability  diagram  shown  in  Figure  6.  The 
boundary  between  dynamic  stability  and  instability  is  given  by  the 
equation: 


I 

-  =  Sd  (2  -  Sd). 
S3 


(3.5) 


This  equation  can  be  used  in  conjunction  with  equation  (3.3)  to  deter¬ 
mine  the  stability  bounds  on  the  projectile  spin  rate  (P).  The  follow¬ 
ing  conclusions  can  be  made  by  referring  to  Figure  6: 


(1)  A  dynamically  stable  projectile  must  be  gyroscopically 

1 

stable  ( i .e. ,  —  <  1 ) ; 

S9 

(2)  If  $d  lies  in  the  interval  (0-2),  a  statically  unstable 
projectile  can  be  stabilized  by  spinning  it  at  a  suffi¬ 
ciently  high  rate  and  a  statically  stable  projectile  is 
always  dynamically  stable; 
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(3)  If  S^  lies  outside  this  interval,  a  statically  unstable 
projectile  cannot  be  spin-stabilized.  In  fact,  a  statically 
stable  projectile  can  be  made  dynamically  unstable  by  spin. 

It  should  be  pointed  out  that  the  results  shown  in  Figure  6  assume  that 
S^  does  not  depend  on  the  spin  rate.  If  there  is  a  strong  dependence 
of  the  magnus  moment  on  spin  then  the  method  of  determining  stability 
must  be  modified  [2]. 

The  effects  of  the  various  coefficients  in  equation  (3.4)  on 
the  dynamic  stability  factor  (S^)  have  been  studied  by  Platou  [SJ.  He 
presents  a  series  of  graphs  that  can  be  used  as  a  guide  to  determine 
the  influence  of  projectile  design  changes  on  the  dynamic  stability. 

In  order  to  apply  the  stability  criterion  described  above,  the 
mass  distribution  (moments  of  inertia,  radii  of  gyration)  and  the  aero¬ 
dynamic  coefficients  in  equations  (3.3)  and  (3.4)  must  be  known.  As 
previously  mentioned,  determination  of  the  aerodynamic  coefficients  for 
self-forging  fragments  by  experimental  wind  tunnel  testing  is  difficult 
and  therefore  numerical  techniques  are  used  for  this  purpose.  A  com¬ 
pendium  of  methods  (DATCOM)  for  determining  aerodynamic  coefficients 
for  subsonic  to  hypersonic  velocity  regimes  has  been  assembled  by  the 
US  Air  Force  [3].  For  the  case  of  hypersonic  velocities  (applicable  to 
self-forging  fragments),  Newtonian  impact  theory  is  used  to  determine 
the  pressure  coefficient  for  any  surface  element.  Fluid  particles  that 
impact  the  surface  are  assumed  to  lose  their  normal  component  of  momen¬ 
tum,  whereas  the  tangential  component  is  preserved.  The  impact  results 
in  a  pressure  coefficient  (C^)  of  the  form 

Cp  =  2  sin''  0  (3.6) 
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where  o  is  the  angle  between  a  tangent  to  the  surface  and  the  direction 
of  the  free  stream  of  fluid  particles.  This  pressure  coefficient  is 
then  used  to  derive  analytical  expressions  for  the  aerodynamic  coeffi¬ 
cients  for  a  particular  shape  of  projectile.  The  DATCOM  method  uses 
design  charts  and  empirical  relationships  to  determine  the  aerodynamic 
coefficients  at  hypersonic  velocities  for  projectile  shapes  composed  of 
one  or  more  cone  frustums  with  or  without  a  spherical  nose.  As  will  be 
shown  later  these  shapes  can  be  used  for  self-forging  fragments. 


The  method  will  now  be  illustrated  for  the  derivative  of  the 

normal -force  coefficient  (C^  )  which  is  given  by 

a 


(3.7) 


where  dfl  and  d^  are  the  base  diameters  of  the  segment  and  projectile, 
respectively.  To  apply  this  equation,  the  body  is  divided  into  m 
segments,  the  first  segment  being  either  a  spherical  nose  or  a  cone 
frustum,  and  each  succeeding  segment  being  a  cone  frustum.  The  para¬ 
meter  C jj  for  a  spherical  nose,  based  on  its  base  area,  is  obtained 

from  Figure  7,  and  for  a  cone  frustum,  based  on  the  base  area  of  the 

>  i 

specific  segment,  is  obtained  from  Figure  8.  The  ratio  (dn/db)2  refers 
to  the  base  area  of  the  configuration.  By  evaluating  for  each 

i.t  i 

segment  and  then  applying  equation  (3.7)  the  derivative  of  the  normal - 
force  coefficient  (C..  )  for  the  complete  projectile  is  obtained.  The 

ii 

•  t 

metnod  applied  to  tne  remainder  of  the  aerodynamic  coefficients  (except 
for  the  derivative  of  the  magnus  moment  coefficient  (C^  ))  is  given  in 

p« 

Appendix  A.  The  derivative  of  the  magnus  moment  coefficient  (C  )  can 

V 

be  determined  from  the  following  empirical  equation  [ 5 J : 
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26.3  (l/d)*L 


(3.8) 


where  1  and  d  are  the  projectile  length  and  diameter,  respectively,  L 
is  the  distance  from  the  magnus  center  of  pressure  to  the  center  of 
gravity  and  is  the  Reynolds  Number  based  on  1. 


All  of  the  information  previously  described  in  this  section  can 
be  incorporated  into  a  computer  program  for  calculating  aerodynamic 
stability  of  self-forging  fragments. 


4.0  SFFSTAB:  AERODYNAMIC  STABILITY  PROGRAM 


A  computer  program  SFFSTAB  has  been  developed  for  calculating 
the  aerodynamic  stability  of  a  self-forging  fragment.  In  tnis  program, 
physical  properties  such  as  volume,  mass,  and  moments  of  inertia  are 
first  calculated.  Formulas  for  these  properties  for  hollow  spherical 
segments  and  hollow  truncated  cones  have  been  derived  and  are  given  in 
Appendix  B.  Next,  the  aerodynamic  coefficients  are  computed  usiny  the 
DATCOM  method  (Appendix  A).  Gyroscopic  and  dynamic  stability  are 
calculated  and  the  projectile's  stability  characteri sties  are  deter¬ 
mined  from  the  stability  criteria  illustrated  in  Figure  6.  If  the 
fragment  is  dynamically  unstable,  but  can  be  spin  stabilized,  then  a 
minimum  fragment  spin  rate  for  stabilization  is  calculated.  Since  the 
fragment  is  initially  formed  from  a  saucer-shaped  liner,  it  is  neces¬ 
sary  to  translate  the  fragment  spin  rate  into  a  spin  rate  for  the 
initial  device  configuration.  This  is  simply  done  by  dividing  tne 
liner  into  segments  and  applying  conservation  of  angular  momentum  (see 
Figure  9). 

As  an  example  of  the  application  of  the  SFFSTAB  computer  pro¬ 
gram  consider  the  self -forging  fragment  radiograph  shown  in  Figure  10. 
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As  illustrated  a  five  segment  model  was  used  to  represent  the  fragment. 
The  objective  of  the  calculation  was  to  determine  the  stability  char¬ 
acteristics  of  tin’s  fragment  and  to  determine  if  it  could  be  spin  stab¬ 
ilized.  The  actual  output  from  the  SFFSTAU  program  is  given  as  the 
sample  problem  in  Appendix  C.  For  the  solid  fragment  shape  shown  in 
Figure  10,  the  calculated  minimum  spin  rate  for  stabilization  was 
ob2B0  RPM.  Based  on  conservation  of  angular  momentum  this  corresponds 
to  an  initial  SFF  charge  spin  rate  of  BB60  RPm. 

SFFSTAB  is  a  useful  tool  for  conducting  parametric  studies  for 
different  fragment  shapes.  The  dependence  of  stability  on  the  exterior 
shape  is  illustrated  in  Figure  11,  where  the  fragment  spin  rate 
required  for  dynamic  stability  as  a  function  of  the  flare  angle  0  is 
shown  for  typical  fragment  shapes.  The  interior  profile  also  has  a 
significant  effect  on  the  dynamic  stability.  Some  typical  results 
obtained  for  fragments  with  similar  exterior  profiles  but  varying 
degrees  of  hollowness  are  shown  in  Figure  12.  Notice  that  as  the  frag¬ 
ment  shape  is  changed  from  hollow  to  completely  solid  the  spin  rate 
required  for  dynamic  stability  increases. 


5.0  SUMMARY 

In  summary,  an  aerodynamic  stability  criterion  which  combined 
gyroscopic  and  dynamic  stability  has  been  presented.  A  technique  for 
calculating  the  aerodynamic  coefficients  was  used  together  with  the 
stability  criterion  for  calculating  the  stability  of  self-forging  frag¬ 
ments.  A  computer  program  (SFFSTAB)  was  developed  using  this  informa¬ 
tion  and  sample  problems  illustrating  the  effects  of  parameters  such  as 
interior  and  exterior  fragment  profiles  have  been  presented.  This 
report  provides  the  complete  documentation  required  to  use  the  program 
SFFSTAB. 
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RADIOGRAPHS  OF  UNSTABLE  FRAGMENT 
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Figure  3 

TECHNIQUES  FOR  SELF- FORGING  FRAGMENT  STABILIZATION 
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Figure  4 

AERODYNAMIC  FORCES  ON  A  BODY 
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Figure  9 

APPROXIMATION  OF  LINER  GEOMETRY  BY  SEGMENTS 
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EXPERIMENT  01 
RADIOGRAPH 


AERODYNAMIC 
MODEL  FOR 
SFFSTAB 


REQUIRED  FRAGMENT  SPIN  RATE  =  6940  (RAD/SEC)  =  66280  (  RPM ) 
REQUIRED  CHARGE  (SFF)  SPIN  RATE  =  404  (RAD/SEC)  =  3860  (RPM) 


Figure  10 

SFFSTAB"  RESULTS  FOR  EXPERIMENTALLY 
OBSERVED  FRAGMENT 
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FRAGMENT  SPIN  RATE  ( RPM  x  10000) 


Figure  11 

FRAGMENT  SPIN  RATE  REQUIRED  FOR 
DYNAMIC  STABILITY  VS  FLARE  ANGLE  0 
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FRAGMENTS  WITH  DIFFERENT  INTERIOR  PROFILES 
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APPENDIX  A  -  DATCOM  METHOD 
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DATCUM  METHOD  FUR  HYPERSONIC  FLOW 


The  DATCUM  Method  [3]  uses  the  results  of  Newtonian  impact 
theory.  Design  charts  and  equations  are  presented  in  Reference  3,  for 
determining  the  aerodynamic  characteristics  at  hypersonic  speeds  of 
projectile  shapes  composed  of  one  or  more  cone  frustums  with  or  without 
a  spherical  nose.  The  coefficients  required  for  gyroscopic  and  dynamic 
stability  equations  are  given  in  this  Appendix.  The  relevant  coeffi¬ 
cients  are: 


1) 

CN  : 

NORMAL-FORCE-CURVE  SLOPE 

ii) 

a 

Cm  : 

PITCHING  MOMENT  CURVE  SLOPE 

iii) 

a 

CD  : 

DRAG  COEFFICIENT 

iv) 

V 

PITCHING  DERIVATIVE 

v) 

q 

Sn  : 

PITCHING  DERIVATIVE 

vi ) 

q 

V* 

300Y  ACCELERATION  DERIVATIVE 

a 


The  information  contained  in  this  Appendix  has  been  extracted  from 
Reference  3. 
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i)  CN  :  NORMAL-FORCE-CURVE  SLOPE 
at 


The  normal -force-curve  slope  of  a  body  composed  of  one  or  more 
cone  frustums  with  or  without  a  spherical  nose,  based  on  tne  body  base 
area,  is  given  by 


dC 


N 


da 


m 

i.  (CN  >n 
n=l  a 


/dn\ 


IV 


( A.  1 ) 


To  apply  this  equation  the  body  is  divided  into  m  segments,  the 

first  segment  being  either  a  spherical  nose  or  a  cone  frustum,  and  each 

succeeding  segment  a  cone  frustum.  The  normal -force-curve  slope  of  a 

spherical  nose,  based  on  area,  is  obtained  from  Figure  A-l.  The 

normal -force-curve  slope  of  a  cone  frustum,  based  on  the  base  area  of 

the  specific  segment,  is  obtained  from  Figure  A-2.  (Note  that  a 

cylinder  is  considered  a  cone  frustum  with  Q  =  0  and  a/d  =  1.0,  and 

that  C^  =  0  by  Newtonian  impact  theory.)  The  ratio  ( dn /db ) 2  refers 
a 

the  normal -force-curve  slope  to  the  base  area  of  the  configuration. 
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Figure  A-1 

NORMAL  FORCE -CURVE  SLOPE  FOR  SPHERICAL  SEGMENTS 
(Reproduced  from  Figure  4.2.1. 1 -23  of  Reference  [3]) 
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i 1 )  CM  :  PITCHING  MOMENT  CURVE  SLOPE 
a 


The  procedure  for  computing  the  total  pi tchi ng-moment-curve 
slope  for  a  complex  body  is  given  in  the  following  steps.  The  moment 
values  for  each  individual  segment  of  a  multiple  cone-frustum  body  with 
or  without  a  spherical  nose  are  referred  to  a  moment  axis  at  the  front 
face  of  that  particular  segment,  and  are  based  on  the  product  of  the 
base  area  and  base  diameter  of  that  particular  segment. 

Step  1:  Compute  C^  for  each  body  segment  about  its  own  front 

a 

face,  using  Figures  A-3  and  A-4. 

Step  2:  Transfer  the  individual  moment  slopes  to  a  common 

reference  axis  by  applying  the  following  moment  trans¬ 
fer  equation  to  each  body  segment. 

n 

Cm  =  cni  +  -  cn  (Per  radian)  (A. 2) 

x  a  d  a 

where 

C^  is  the  normal -force-curve  slope  of  the  individual 

a 

cone-frustum  or  spherical  nose  segment,  based  on 
its  own  base  area,  from  Figures  A-2  and  A-l, 
respectively. 

d  is  the  base  diameter  of  the  individual  cone- 
frustum  or  spherical  nose  segment. 

n  is  the  distance  from  the  front  face  of  a  given 
segment  to  the  desired  moment  reference  axis  of 
the  configuration,  positive  aft. 
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<1 


is  the  pitching-moment-curve  slope  of  an  indi 


vidual  segment  from  Figure  A-3  for  cone  frustums 
and  from  Figure  A-4  for  spherical  nose  segments. 


C‘  is  based  on  the  product  of  the  base  area  and 
a 


the  base  diameter  of  the  individual  segment. 


Cm  is  the  pitching-moment-curve  slope  of  an  indi - 

a 


vidual  segment  based  on  the  product  of  tne  base 
area  and  base  diameter  of  the  individual  segment 
and  referred  to  a  common  reference  axis. 


Step  3:  The  transferred  pitching-moment-curve  slopes  of  the 
individual  body  segments  are  then  converted  to  a 
common  basis  by 


Cm  =  l  (Cm  »n  “  (Per  radl'an) 
a  n=!  a  db  I 


(A. 3) 


where  the  subscript  n  refers  to  an  individual  segment 


of  m  segments,  and  C  is  referred  to  a  common  refer- 

a 


ence  axis  and  is  based  on  the  product  of  the  area  and 
diameter  of  the  base  of  the  configuration  Sbdb> 
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SEMICONE  ANGLE,  6  (deg) 


C'm&  FOR  CONE  FRUSTUMS 
(Reproduced  from  Figure  4.2. 2.1 -25a  of  Reference  [  3  ] ) 
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NOSE  BLUNTNESS,  2<’/ds 


Figure  A- 4 

C'm  FOR  SPHERICAL  SEGMENTS 

m0 

(Reproduced  from  Figure  4. 2.2.1 -25b  of  Reference  [3D 
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iii)  CQ :  DRAG  COEFFICIENT 


The  zero-lift  drag  (based  on  the  maximum  frontal  area)  of  a 
body  composed  of  one  or  more  cone  frustums  with  or  without  a  spherical 
nose  is  estimated  by  adding  the  pressure-drag  coefficient  of  each  seg¬ 
ment  to  the  body  skin-friction  drag  coefficient. 


(A. 4) 


The  procedure  to  be  followed  in  evaluating  this  equation  is  as 
follows: 

Step  1:  Divide  the  body  into  m  segments,  the  first  segment 

being  either  a  spherical  nose  or  a  cone  frustum,  and 

each  succeeding  segment  a  cone  frustum.  The  pressure- 

drag  coefficient  for  a  spherical  nose  or  cone  frustrum 

is  obtained  from  Figures  A-5  and  A-6,  respectively. 

The  pressure-drag  coefficients  for  the  remainder  of 

the  segments  are  obtained  from  Figure  A-7.  The 

pressure-drag  coefficients  are  based  on  the  base  area 

of  the  specific  segment.  The  ratio  (d  /d  )2  refers 

n  max 

the  pressure-drag  coefficients  to  the  maximum  body 
frontal  area. 

Step  2:  Obtain  the  body  skin- friction  drag  coefficient  by 

Cf  $ 
c  S 

C  =  1.02  Cf -  (A. 5) 

f  inc  Cf  Sg 
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where 

is  the  incompressible  (M  =  0)  turbulent,  flat- 

plate  skin-friction  coefficient,  including  rough¬ 
ness  effects,  as  a  function  of  Reynolds  number 
based  on  the  total  length  of  the  body  Jig.  This 
value  is  obtained  from  Figure  A-8. 

i  is  the  reference  length  in  inches. 

k  is  the  surface-roughness  height  in  inches;  it 
depends  upon  surface  finish.  Representative 
values  for  this  parameter  can  be  obtained  from 
Table  A-l. 

The  ratio  1/ k  is  computed  and  Figure  A-9  is  used  to 
obtain  the  cutoff  Reynolds  number.  If  the  cutoff 
Reynolds  number  is  greater  than  the  computed  Reynolds 
number  for  the  specific  configuration,  the  value  of 
is  obtained  from  Figure  A-8  at  the  computed  Reynolds 
number.  If  the  cutoff  Reynolds  number  is  less  than 
the  computed  Reynolds  number,  the  value  of  is 
obtained  from  Figure  A-8  at  the  cutoff  Reynolds 
number. 

is  the  ratio  of  compressible  to  incompressible 

skin-friction  coefficient  obtained  from 
Figure  A-10. 
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—  is  the  ratio  of  the  bocjy  wetted  area  to  maximum 
SB 

body  frontal  area. 

If  this  method  is  applied  at  Mach  numbers  low  enough 
so  that  the  base  drag  is  significant,  the  base  drag 
should  be  added  to  the  results  obtained.  Unfortun¬ 
ately,  the  only  base-drag  coefficient  results  avail¬ 
able  which  are  compatible  with  the  Newtonian-theory 
results  (restricted  to  bodies  with  forward  facing 
slopes  or  cylinders,  in  which  case  the  Newtonian 
results  are  equal  to  zero)  are  those  for  cylindrical 
afterbodies.  The  pressure-drag  coefficient  for  cylin¬ 
drical  afterbodies  is  presented  in  Figure  A-ll  for 
M  <  10. 
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Figure  A-7 

DRAG  FORCE  COEFFICIENT  DUE  ONLY  TO  THE  INCLINED  SIDES 
OF  A  CONE  FRUSTUM  CALCULATED  BY  NEWTONIAN  THEORY. 
CD  IS  BASED  ON  BODY  BASE  AREA  Sb- 
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TURBULENT  MEAN  SKIN  -  FRICTION  COEFFICIENT  ON  AN  INSULATED  FLAT  PLATE 
(Reproduced  from  Figure 4.1. 5.1 -26  of  Reference  [3]) 
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Table  A-1 

REPRESENTATIVE  VALUES  OF  SURFACE  ROUGHNESS  HEIGHT 
(Reproduced  from  Table  4.1.5.1  -  A  of  Reference  [3]) 


TYPE  OF  SURFACE 

EQUIVALENT  SAND  ROUGHNESS 

k  (inches) 

Aerodynamically  Smooth 

0 

Polished  Metal  or  Wood 

0.02  -  0.08  x  10" 3 

Natural  Sheet  Metal 

0.16  x  10'3 

Smooth  Matte  Paint,  Carefully  Applied 

0.25  x  10  3 

Standard  Camouflage  Paint,  Average  Application 

0.40  x  10' 3 

Camouflage  Paint,  Mass- Production  Spray 

1.20  x  10' 3 

Dip -Galvanized  Metal  Surfaces 

6  x  10  3 

Natural  Surface  of  Cast  Iron 

10  x  10' 3 
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MACH  NUMBER 


2  4  8  1  2  46  1  2  4  6  1  2  46 


CUTOFF  REYNOLDS  NUMBER  (Rf) 


Figure  A -9 

CUTOFF  REYNOLDS  NUMBER 
(Reproduced  from  Figure  4.1.5.1 -27  of  Reference  [3]) 
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(Adapted  from  Figure  4.2.3.1  -  68  of  Reference  [3]) 


UNCLASSIFIED 


SM  1194 


BASE  DRAG  COEFFICIENT  FOR  BODIES  OF  REVOLUTION 

WITH  NO  BOATTAIL 

(Adapted  from  Figure 4.2.3.1  - 60  of  Reference  [3]) 
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iv)  CN  :  PITCHING  DERIVATIVE 


Charts  based  on  simple  Newtonian  theory,  are  presented  for 
determining  of  spherical  segments  and  cone  frustums  at  small  angles 

q 

of  attack. 


The  coefficients  of  these  charts  are  referred  to  the  body  base 
area  and  base  diameter  and  to  a  moment  center  at  the  forward  face  of 

the  segment.  By  proper  use  of  the  data  presented,  the  total  C.,  may  be 

q 

determined  for  bodies  composed  of  multiple  cone  frustums  with  or 
without  spherically  blunted  noses. 


The  Newtonian  value  of  the  derivative  for  a  complex  body  is 

q 

obtained  as  follows: 


Step  1 :  Compute  C^  for  each  bocjy  segment  about  i  ts  front  face 

q 

using  Figures  A-12  and  A-13. 


Step  2:  Transfer  the  individual  derivatives  C'  to  a  common 

N 

x  .  q 

reference  axis  by  appling  the  following  transfer  equa¬ 
tion  to  each  body  segment 


CN  CN  '  2  (“)CN 

q  q  da 


:a.6) 


where 


Cjg  is  the  normal -force-curve  slope  for  each  segment 

a 

based  on  individual  base  areas. 
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CjJj  is  the  pitching  derivative  for  each  segment  based 

q 

on  individual  base  areas  and  base  diameters  and 
referred  to  a  moment  center  at  the  forward  face  of 
the  segment,  from  Figures  A-12  and  A-13. 

n  is  the  distance  from  the  face  of  a  given  frustum 
to  the  desired  moment  reference  axis  of  the  con¬ 
figuration,  positive  aft. 


d  is  the  base  diameter  of  a  given  frustum. 


Step  3:  The  transferred  derivatives  of  the  individual  body 

segments  are  converted  to  a  common  reference  area  and 
diameter  and  added.  The  total  derivative  is  given  by 


1  (c»  >„ 

n=l  q 


(A.7) 
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CN  BASED  ON  BASE  AREA 
q  AND  BASE  DIAMETER 

T  I - 1 - I - I - 1 - 1-  . 1 - 1 - — 1 

1  .2  3  .4  5  .6  7  8  .9  1.0 

NOSE  BLUNTNESS  121’ /ds) 


Figure  A  -13 

PITCHING  DERIVATIVE  C'N  FOR  SPHERICAL  SEGMENTS 
(Reproduced  from  Figure  7.2. 1.1  9b  of  Reference  131) 
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v)  C  :  PITCHING  DERIVATIVE 

m 


Charts  based  on  simple  Newtonian  theory  are  presented  for 
determining  C^  of  spherical  segments  and  cone  frustums  at  small  angles 

q 

of  attack.  The  coefficients  of  these  charts  are  referred  to  the  base 
area  and  the  square  of  the  base  diameter  and  to  a  moment  center  at  the 
forward  face  of  the  segment.  By  proper  use  of  the  data  presented,  the 

total  C^  may  be  determined  for  bodies  composed  of  multiple  cone 

q 

frustums  with  or  without  spherically  blunted  noses. 

The  Newtonian  value  of  the  stability  derivative  C  for  a 

mq 

complex  body  is  obtained  as  follows: 


Step  1:  Compute  C  '  for  each  body  segment  about  its  own  front 

q 

face  using  Figures  A-14  and  A-16  if  the  body  has  a 
spherically  blunted  nose,  and  Figures  A-14  and  A- 15  if 
the  body  nose  is  a  cone  frustum. 


Step  2:  Transfer  the  individual  derivatives  C'  to  a  common 

fil 

q 

moment  center  by  applying  the  following  axis  transfer 
equation  to  each  body  segment: 


n 


'm 


=  C 


„  -  2  -  C  ‘ 

m  ,  m 

q  da 


n 

+  -  C 


d  Nq 


-  2(-|  C 
d  , 


M 


(A. 8) 


where 


CfJ  is  the  normal  -  force-curve  slope  for  each  segment 
1 

based  on  individual  base  areas. 
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CjJj  is  the  pitching  derivative  for  each  segment  based 
d 

on  individual  base  areas  and  base  diameters  and 
referred  to  a  moment  center  at  the  forward  face  of 
the  segment. 

is  the  pitching-moment-curve  slope  for  each  seg- 

a 

ment  based  on  individual  base  areas  and  base  dia¬ 
meters  and  referred  to  a  moment  center  at  the 
forward  face  of  the  segment. 


CjJ,  is  the  pitching  derivative  for  each  body  segment 

q 

based  on  individual  base  areas  and  the  square  of 
base  diameters  and  referred  to  a  moment  center  at 
the  forward  face  of  the  segment.  If  a  complex 
body  consists  of  combinations  of  cone  frustums, 
the  derivative  for  the  first  frustum  must  be 
obtained  from  Figure  A-15,  which  accounts  for  the 
front  face  being  exposed  to  the  air  stream.  If 
the  body  has  a  spherically  blunted  nose,  the 
derivative  of  the  nose  is  obtained  from 
Figure  A-16.  For  subsequent  frustums  the  deriva¬ 
tives  are  obtained  from  Figure  A-14. 


Step  3:  The  transferred  derivatives  of  the  individual  body 
segments  are  added  after  being  converted  to  a  common 
reference  area  and  squared  diameter.  The  total  deriv¬ 
ative  of  the  individual  body  segments  is  given  by 


m 


(A. 9) 
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The  boc(y  contribution  to  the  derivative  Cm#  in  the  hypersonic 

a 


speed  range  is  equal  to  zero  when  determined  by  the  Newtonian  Impact 
Theory. 
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j  APPENDIX  B  -  SEGMENT  PROPERTIES 

l 

i 

1 


I  UNCLASSIFIED 

i 


UNCLASSIFIED 


PROPERTIES  OF  A  SPHERICAL  SEGMENT 


Volume 


R3  2r3 

V  =  n  [Rr2  +  hR2  -  hr2  - 

3  3 


Centroid 


3  [  (R2-r2)(4hR-2h2)  -  (R2-r2 )2 , 

y  = - 

4  [3h  {R2-r2 )  +  3Rr2  -  R3  -  2r3j 


Axial  Moment  of  Inertia  {Taken  about  Y-axis) 


3  2 

M  [—  (R b -r 5 )  +  (R-h){r4-R‘i)  +-  (R-h) 3  (R2-r2 )  j 
15  3 


M  R3  2r3 

2  iRr2  +  hR2  -  hr2 - ] 

3  3 

Where  M  =  Mass  =  Density  x  Volume 
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<  I 


UNCLASSIFIED 


SM  1194  B-4 


HOLLOW  TRUNCATED  CONE 


r  i  =  INNER  RADIUS  AT  SMALL  END  OF  TRUNCATED  CONE 

r2  =  OUTER  RADIUS  AT  SMALL  END  OF  TRUNCATED  CONE 

R !  =  INNER  RADIUS  AT  LARGE  END  OF  TRUNCATED  CONE 

R  2  =  OUTER  RADIUS  AT  LARGE  END  OF  TRUNCATED  CONE 

h  =  SEGMENT  HEIGHT 
y  =  CENTROID 


,*V  v 


mmm 
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PROPERTIES  OF  A  HOLLOW  TRUNCATED  CONE 


Volume 


rch 

v  =  _  [R/  -  Rj2  +  r22  -  rt2  +  R2r2  -  R 1  r ^ 


Centroid 


h  LR22  -  Rt2  +  3r22  -  3rx2  +  2R2r2  -  ZR^J 
4  [R22  -  Rt2  +  r22  -  rt2  +  R2r2  -  R^J 


Axial  Moment  of  Inertia  (Taken  about  Y-axis) 


3M 


I A  =  —  [  ( R2  **  +  R23r2  +  R22r/  +  R2r23  +  r24^ 


10 


-  (Rj-4  +  R13r1  +  R^r/  +  R^3  +  r^)] 
/  lrz2  .  Rj2  +  .  r ^ 2  +  R2r2  -  R^J 


Where  M  =  Mass  =  Density  x  Volume 
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s 


ti 


■S- 


& 


& 
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TRANSVERSE  MOMENT  OF  INERTIA 


(Taken  Parallel  to  X-axis  and  Through  the  Centroid) 


3M  (R2-r2)“ 

IT  =  l—  {R24  "  2R23  (R2"r2>  +  2R22  (R2 ~*^2 ^ 2  -  + - 


(Ri"r 1 )4 

-  Rj4  +  2R j 3  (Rj-rj)  -  2Rj2  (R1-r1)2  +  R1(Rl-r1)3 - ) 


3M  (h2  { 


R22  (R2'r2^2  (R2-r2)  Rx2  (R^j)2  (Rj-rJ 


-  R 


+  R, 


'2 - - *  *  I  - 

3  5  2  3  5  2 


(R2-r2)2  (R,-r i  )2 

+  y2  {R22  + - R2  (R2-r2)  -  Rt2 - +  RL  (R1-r1 ) } 


^R2-r2^2  (R2"^2)  ^Rl“rJ2  (Rl'rJ 


+  yh  {-R22 - - — - —  +  4R2  — - — —  +  R^  + 


4R , 


2  ‘  3 


!  lR22  -  V  ♦  <■/  -  r*  *  R2r2  -  R,rJ 


V/ 
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APPENDIX  C 


-  FlQy  CHART  AND  FORTRAN  LISTING  OF  SFFSTAb 
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SFFSTAB:  FLOWCHART 


I 


READ  AND  WRITE  INPUT  DATA 


CALCULATE  PROPERTIES  FOR  EACH  SEGMENT 


W  COEFFI 


W  COEFF3 


U  COEFF4 


Cu  NORMAL-FORCE  CURVE 

N 

a 

SLOPE  FOR  SPHERICAL  SEGMENTS 


COEFF2 


CN  NORMAL-FORCE  CURVE 

a 


SLOPE  FOR  CONE  FRUSTUMS 


C'  PITCHING  MOMENT 
in 

a 


CURVE  SLOPE  FOR  CONE  FRUSTUMS 


C'  PITCHING  MOMENT 
in 

a 

CURVE  SLOPE  FOR  SPHERICAL  SEGMENTS 


COEFF7  |  Cf  TURBULENT  MEAN  SKIN 


FRICTION  COEFFICIENT  ON  AN 
INSULATED  FLAT  PLATE 


COEFFH  ]~C  BASE  DRAG 


.oefucilnt  with  no  boattail 


k 

COEFFS 

j  Cf(./Cf  COMPRESSIBILITY 

ON  TURBULENT  SKIN  FRICTION 

1 

EFFECT 

I 

J 

4 

-* 

C0EFF6 

R  CUTOFF  REYNOLDS 

4 

NUMBER 

F  I bURE  ^-1  >PPsTAfct  FLOW  LHART 


SUBROUTINES  TO 

GENERATE  CUBIC 
SPLINE  COEFFICIENTS 
FOR  DATCOM  GRAPHS 


JNCLASb IF IED 
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§ 


I 
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FIGURE  C-2  SFFSTAB  FLOWCHART  CONT'O 
UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 
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i  Mt.UW  C-S  iHSTAU  HUM  CHART  CONT'j 

\ 

i 

f 
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IGURt  C-b  SFFSTAB  FLOWCHART  CONT'J 
UNCLASSIFIED 


UNCLASSIFIED 
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C  PROGRAM  TO  CALCULATE  GYROSCOPIC  AND  DYNAMIC  STABILITY 
C 

DIMENSION  ISEG(IOO) ,R2(100) , R1 ( 100) , H( 100 ) ,RR2(100) , RR1 ( 100 ) 
DIMENSION  V(100) .AMASS (100) ,YB(100) ,AI(100) ,TI(100) ,CNAP(10Q) 
DIMENSION  CMAP( 100) .THETA(IOO) ,AD(100) ,DBCS( 100) ,DCBCS( 100) 
DIMENSION  DB(IOO) ,DBFS(100) ,DBFSRA(1Q0) ,CMA(100) ,CMAL(100) 
DIMENSION  COP ( 1 00 ) ,AREA(100) , CNQP ( 100 ) ,CMQP(100) ,CMQ(100) 
C0MM0N/XY1/X1( 11) ,Y1(11) 

C0MM0N/XY2/X2( 11 ,10),Y2(11,10) 

COMMON/XY3/X3( 11 , 10) , Y3( 11 ,10) 

COMMON/X Y4/X4 ( 1 3 ) , Y4 ( 1 3 ) 

COMMON/XY5/X5(8) ,Y5(8) 

C0MM0N/XY6/X6( 5) ,Y6(5) 

COMMON/XY7/X7(26) ,Y7(26) 

C0MM0N/XY8/X8(8),Y8(8) 

COMMON/XYB/XXB ( 6 ) , YYB ( 6 ) 

C0MM0N/BCD1/B1(11),C1(11),D1(11) 

COMMON/BCD2/B2( 11, 10), C2 (11, 10), D2 (11, 10) 

COMMON/BCD3/B3( 11, 10), C3( 11, 10), 03(11,10) 

C0MM0N/BCD4/B4( 13) ,C4( 13) ,D4( 13) 

COMMON/BCD5/B5(8) ,C5(8) ,D5(8) 

C0MM0N/BCD6/B6( 5) ,C6(5) ,06(5) 

COMMON/BCD7/B7(26) ,C7(26) ,D7(26) 

C0MM0N/BCD8/B8(8) ,C8(8) ,D8(8) 

COMMON/BCD3B/B3B(6) ,C3B(6) ,D3B(6) 

CHARACTER*6  TITLE( 15) 

C 

Q* ******  *********************  ****************************** 

C 

C  ITYPE  =  0  -  SPHERICAL  SEGMENT  AND  TRUNCATED  CONES  ARE  USED 
C  1  -  ONLY  TRUNCATED  CONES  ARE  USED 

C  DEN  =  DENSITY  OF  MATERIAL  (MASS/VOLUME) (LB  SEC**2/IN**4) 

C  N:EG  =  NUMBER  OF  SEGMENTS 

C  DBCS  =  DISTANCE  FROM  BASE  OF  BODY  TO  CENTROID  OF  SEGMENT 
C  TVOL  =  TOTAL  VOLUME  OF  B0DY(IN**3) 

C  TMASS  =  TOTAL  MASS  OF  BODY(LB  SEC**2/IN) 

C  CGOB  =  CENTER  OF  GRAVITY  OF  BODY  (RELATIVE  TO  BASE) 

C  AIB  =  AXIAL  (POLAR)  MOMENT  OF  INERTIA  OF  THE  BODY(LB  SEC**2  IN) 

C  TIB  =  TRANSVERSE  MOMENT  OF  INERTIA  OF  THE  BODY(LB  SEC**2  IN) 

C  DCBCS  =  DISTANCE  FROM  CENTROID  OF  BODY  TO  CENTROID  OF  SEGMENT 
C  DB  =  BASE  DIAMETER  OF  THE  SEGMENT 
C  CNAT  =  TOTAL  CNA  FOR  THE  BOOY 

C  CMAT  =  TOTAL  CMA  FOR  THE  BODY  (BASED  ON  AREA  AND  BASE  DIAMETER) 

C  CMATL  =  TOTAL  CMA  FOR  THE  BODY  (BASED  ON  AREA  AND  BODY  LENGTH) 

C  RA  =  DISTANCE  FROM  MOMENT  REFERENCE  AXIS  TO  BASE  OF  BODY 
C  ( IF  RA  =  0.0  THE  C  OF  G  OF  THE  BODY  IS  USED) 

C  DBFS  =  DISTANCE  FROM  BASE  OF  BODY  TO  FRONT  FACE  OF  SEGMENT 

C  DBFSRA  =  DISTANCE  FROM  FRONT  FACE  OF  SEGMENT  TO  MOMENT  REFERENCE 

C  AXIS 

C  CMA  =  CMA  FOR  A  SEGMENT  ,  REFERRED  TO  A  COMMON  REFERENCE  AXIS 
C  (BASED  ON  AREA  AND  BASE  DIAMETER) 

C  CMAL  =  CMA  FOR  A  SEGMENT  ,  REFERRED  TO  A  COMMON  REFERENCE  AXIS 
C  (BASED  ON  AREA  ANO  SEGMENT  LENGTH) 

C  BL  =  LENGTH  OF  BOOY 


UNCLASSIFIED 


UNCLASSIFIED 


C-8 


C  XCP  =  CENTER  OF  PRESSURE  AS  A  FRACTION  OF  BODY  LENGTH 
C  XH  =  MOMENT  CENTER  LOCATION  AS  A  FUNCTION  OF  BODY  LENGTH 
C  RHO  =  AIR  DENSITY(LB  SEC**2/IN**4) 

C  W  =  SPIN  RATE  (RAD/SEC) 

C  VEL  =  PROJECTILE  VELOCITY( IN/SEC) 

C  ISG  =  1  -  CALCULATE  W 

C  0  -  CALCULATE  SG 

C  SG  =  GYROSCOPIC  STABILITY  FACTOR 

C  CDP  =  NEWTONIAN  DRAG  COEFFICIENT 

C  BLUNT  =  NOSE  BLUNTNESS  FOR  A  SPHERICAL  SEGMENT 

C  ESR  =  EQUIVALENT  SAND  ROUGHNESS(IN) 

C  RE  =  REYNOLDS  NUMBER  (BASED  ON  L  ) 

C  RES  =  SAME  AS  RE 
C  CRE  =  CUTOFF  REYNOLDS  NUMBER 
C  AR  =  ADMISSIBLE  ROUGHNESS 
C  U  =  ABSOLUTE  VISCOSITY(LB  SEC/IN**2) 

C  CF  =  INCOMPRESSIBLE  FLAT  PLATE  SKIN  FRICTION  COEFFICIENT 
C  CFCCF  =  RATIO  OF  COMPRESSIBLE  TO  INCOMPRESSIBLE  SKIN  FRICTION 
C  COEFFICIENT 

C  SB  =  BODY  MAXIMUM  FRONTAL  AREA(IN**2) 

C  SS  =  WETTED  AREA  OR  SURFACE  AREA  OF  THE  BODY  EXCLUDING  THE  BASE 

C  CDF  =  BODY  SKIN  FRICTION  DRAG  COEFFICIENT 

C  CDO  =  ZERO  LIFT  DRAG  COEFFICIENT 

C  CDB  =  BASE  DRAG  COEFFICIENT 

C  AREA  =  SUFACE  AREA  OF  A  SEGMENT 

C  CNQP  =  PITCHING  DERIVATIVE  FOR  A  SEGMENT 

C  CMAP  =  PITCHING  MOMENT  CURVE  SLOPE  FOR  A  SEGMENT 

C  CNAP  =  NORMAL  FORCE  CURVE  SLOPE  FOR  A  SEGMENT 

C  CMQP  =  PITCHING  DERIVATIVE  FOR  SEGMENT 

C  CMQ  =  PITCHING  OERVATIVE  TRANSFERRED  TO  A  COMMON  AXIS 

C  CMQT  =  PITCHING  DERIVATIVE  FOR  BODY 

C  CMPA  =  DERIVATIVE  OF  THE  MAGNUS  MOMENT  COEFFICIENT 

C  CMAD  =  ACCELERATION  DERIVATIVE 

C  SO  =  DYNAMIC  STABILITY  FACTOR 

C  AKB  =  DIMENSIONLESS  AXIAL  RADIUS  OF  GYRATION 

C  TKB  =  DIMENSIONLESS  TRANSVERSE  RADIUS  OF  GYRATION 

C  DMCP  =  DISTANCE  FROM  MAGNUS  CENTER  OF  PRESSURE  TO  CENTER  OF  GRAVITY 

C  ISPIN  =  1  -  READ  LINER  DATA 

C  TAM  =  TOTAL  ANGULAR  MOMENTUM  (LB  SEC  IN) 

C  WL  =  SPIN  RATE  OF  LINER  (RAD/SEC) 

C  RPML  =  SPIN  RATE  OF  LINER  (RPM) 

C  NLS  =  NUMBER  OF  LINER  SEGMENTS 
C  TLS  =  THICKNESS  OF  LINER  SEGMENT 
C  HLS  =  HEIGHT  OF  LINER  SEGMENT 

C  RCS  =  RADIUS  TO  CENTROIO  OF  CROSS-SECTIONAL  AREA  OF  LINER  SEGMENT 
C  AMASSL  =  MASS  OF  LINER  SEGMENT 
C  TMASSL  =  TOTAL  MASS  OF  THE  LINER 
C 

£********************************************************** 

C 

C  REAO  AND  WRITE  INPUT  DATA 
C 

READ( 5,900) (TITLE( I), 1=1, 15) 

900  FORMAT( 13A6 ,2A1 ) 
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WRITE ( 6 , 801 )(TITLE( I) ,1=1,15) 

801  F0RMAT(1H  ,///,5X ,13A6,2A1 ,//) 

READ(5,10)DEN 
10  FORMAT(F10.5) 

READ( 5, 15) ITYPE ,NSEG 
15  FORMAT (215) 

IF( ITYPE .EQ. 1)G0  TO  20 

READ( 5 ,25) ISEG( 1) ,R2(1),R1(1),H(1) 

25  FORMAT ( 1 5 , 5X , 3F10 . 5) 

DO  30  1=2 ,NSEG 

READ(5 ,35) ISEG( I ) ,RR2( I) ,RR1(I),R2(I),R1(I),H(I) 

35  FORMAT( 15 , 5X , 5F10 . 5) 

30  CONTINUE 
GO  TO  40 

20  DO  45  1=1 ,NSEG 

READ( 5,35) ISEG( I),RR2(I),RR1(I),R2(I),R1(I),H(I) 

45  CONTINUE 

40  WRITE(6 ,46) 

46  FORMAT(1H1,10X, '*****  GYROSCOPIC  STABILITY  DATA  *****',//) 
WRITE(6,50)DEN, ITYPE, NSEG 

50  F0RMAT(1H  ,5X, 'DENSITY  =' ,F15 .8 , 5X , 1 ITYPE  = ' , 15 ,5X , 

1  'NUMBER  OF  SEGMENTS  =' ,15) 

WR ITE (6,55) 

55  F0RMAT(1H  ,7X, ' I SEG ' ,9X, ' RR2 1 ,12X, ' RR1 ' ,13X, ' R2 ' ,13X, 

1  ' R1 1 , 13X , ' H ' ,/) 

IF( ITYPE .EQ. 1 )G0  TO  60 

WR ITE(6 ,65) ISEG( 1 ) ,R2(1),R1(1),H(1) 

65  F0RMAT(1H  ,5X,I5,30X,3(5X,F10.5)) 

DO  70  1=2, NSEG 

WRITE(6,75)ISEG(I) ,RR2(I) ,RR1(I),R2(I),R1(I),H(I) 

75  F0RMAT(1H  ,5X,I5,5(5X,F10.5)) 

70  CONTINUE 
GO  TO  80 

60  DO  90  1=1, NSEG 

WRITE(6 ,75) ISEG( I ) ,RR2( I),RR1(I),R2(I),R1(I),H(I) 

90  CONTINUE 
80  READ(5,210)RA 
210  F0RMAT(F10.5) 

READ(5,295) ISG 
295  FORMAT( 15) 

WRITE(6,400)RA, ISG 

400  F0RMAT(1H  , 5X ,' REFERENCE  AXIS  =' , F 10 . 5 , 5X , ' ISG  = '  ,15) 

IF( ISG. EQ. 1 )G0  TO  300 
READ(5,305)W,RH0,VEL 
305  FORMAT(3F10.5) 

WRITE(6 ,405)W,RHO,VEL 

405  FORMAT( 1H  ,5X,'W  = ' ,F20 . 10 , 5X , ' RHO  =' ,F20. 15, 5X, 'VELOCITY 
1  F15.5) 

GO  TO  420 

300  READ ( 5 , 305)SG ,RHO, VEL 
WRITE(6,410)SG,RHO,VEl 

410  FORMAT( 1H  ,5X,'SG  = ' ,F10 . 5 , 5X , ' RHO  =' ,F20 . 15 , 5X ,' VELOCITY 
1  F15.5) 
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C  CALCULATE  PROPERTIES  FOR  EACH  SEGMENT 
C 

420  CALL  COEFF1 
CALL  COEFF2 
CALL  COEFF3 
CALL  COEFF4 
CALL  COEFF5 
CALL  COEFF6 
CALL  COEFF7 
CALL  COEFF8 
DO  100  1=1 , NSEG 
IF(I.GT.l)GO  TO  105 
IF(ITYPE.EQ.l)GO  TO  105 

CALL  IHSS(ISEG(1) ,R2(1),R1(1),H(1),DEN,V(1) ,AMASS(1) ,YB(1) , 
1  B , A I ( 1 ) , TI ( 1 ) ) 

YB(1)=B 

BLUNT=H(1)/R2(1) 

CALL  F421123(BLUNT,CNAP(1)) 

CALL  F422125B(BLUNT,CMAP(1)) 

WRITE(6,110)BLUNT,CNAP( 1) ,CMAP(1) 

110  F0RMAT(1H  ,5X,'NOSE  BLUNTNESS  =  ' .F10.5.5X, 'CNAP  =',F10.5, 

1  5X , 1 CMAP  =' ,F10 . 5 ,//) 

GO  TO  100 

105  CALL  IHTC( ISEG( I ) ,RR2( I) ,RR1( I),R2(I),R1(I),H(I),DEN,V(I). 

1  AMASS( I) ,YB(I),AI(I),TI(I)) 

AD( I )=R2( I )/RR2( I ) 

THETA( I )=ATAN( ( RR2( I )— R2( I ) )/H( I ) )/3 . 14159265*180 . 

CALL  F421126(THETA( I ) ,CNAP( I ) ,AD( I ) ) 

CALL  F422125A( THETA( I ) , CMAP( I ) , AD( I ) ) 

WRITE(6,115)THETA( I) ,AD( I ) ,CNAP( I) ,CMAP( I ) 

115  FORMAT( 1H  ,5X,'THETA  =' ,F10.5,5X, 1  A/D  =',F10.5,5X, 

1  'CNAP  =' .F10.5.5X, 'CMAP  =',F10.5,//) 

100  CONTINUE 

CALCULATE  DISTANCE  FROM  BASE  OF  BODY  TO  CENTROID  OF  SEGMENT 

SUM1=H(NSEG) 

NSEGM1=NSEG-1 
DBCS( NSEG )=YB( NSEG) 

DO  120  1=1 ,NSEGM1 
DBCS(NSEG-I)=YB(NSEG-I)+SUM1 
SUM1=SUM1+H(NSEG-I) 

120  CONTINUE 

CALCULATE  THE  TOTAL  VOLUME  AND  TOTAL  MASS  OF  BODY 

TVOL=0.0 
TMASS=0 . 0 
DO  125  1=1, NSEG 
TVOL=TVOL+V( I ) 

TMASS=TMASS+AMASS( I ) 

125  CONTINUE 

C  CALCULATE  CENTER  OF  GRAVITY  OF  BODY  (RELATIVE  TO  BASE) 
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SUM2=0 . 0 
DO  130  1=1 ,NSEG 
SUM2=SUM2+DBCS( I )*V( I ) 

130  CONTINUE 

CG0B=SUM2/TV0L 

C 

C  CALCULATE  AXIAL  (POLAR)  MOMENT  OF  INERTIA 
C 

AIB=0 . 0 

DO  135  1=1 ,NSEG 
A I  B=A  I  B-*-A  I  ( I ) 

135  CONTINUE 
C 

C  CALCULATE  DISTANCE  FROM  CG  OF  BODY  CG  OF  SEGMENT 
C 

DO  140  1  =  1  ,NSEG 

DCBCS( I )=ABS(DBCS( I )-CGOB) 

140  CONTINUE 
C 

C  CALCULATE  TRANSVERSE  MOMENT  OF  INERTIA  USING  PARALLEL  AXIS  THEOREM 
C 

TIB=0 . 0 

DO  145  1=1 ,NSEG 

TIB=TIB+TI(I)>AMASS(I)*DCBCS(I)**2. 

145  CONTINUE 

C 

C  CALCULATE  BASE  DIAMETER  OF  EACH  SEGMENT 
C 

DO  195  1=1 ,NSEG 

I F ( I . GT . 1 ) GO  TO  190 

IF( ITYPE . EQ  1 )G0  TO  190 

DB(1)=2.*(H(1)*(2*R2(1)-H(1)))*V5 

GO  TO  195 

190  DB(I)=2.*RR2(I) 

195  CONTINUE 
WRITE ( 6 , 150 ) 

150  FORMAT ( 1H  ,6X , ' I SEG ' ,8X , 'DECS' . 1 1 X . ' DCBCS ‘ . 12X. 'OB' 

DO  155  1=1 ,NSEG 

WRITE(  6 , 160) ISEG( I) ,DBC$( I ) ,DCBCS( I ) ,DB( P 
160  FORMAT ( 1H  , 5X , 1 5 , 3( 5X , F 10  5 ) ) 

155  CONTINUE 

WR ITE ( 6 , 165 ) TVOL 

165  FORMAT ( 1H  ,5X, 'TOTAL  VOLUME  OF  BODY  ='.F10  51 
WR I TE ( 6 , 170) TMA  S  S 

170  FORMAT ( 1H  ,5X,'TOTAL  MASS  OF  BODY  ='.F10  5) 

WR ITE( 6 , 175) CGOB 

175  FORMAT ( 1H  ,5X, 'CENTER  OF  GRAV I  TV  OF  BODy  PpA’PE  ’•  BA^f 
1  F 1 0 . 5 ) 

WR  I  TE  (  6 , 1 80 )  A I B 

180  F0RMAT(1H  .  5X  .  'AXIAL  MOMENT  OF  INERTIA  V  BOD'*  .FI'  5 


lii 


WRITE ( 6 , 185  )TIB 

185  FORMAT ( 1H  ,SX. 'TRANSVERSE  MQME  N T  INEC'IA  >  BOD> 


1  INCLASS  IF  IEP 


TSBSr-' 


UNCLASSIFIED 


C-12 


CALCULATE  CNA  FOR  THE  BODY 

CNAT=0 . 0 
DO  200  1=1 .NSEG 

CNAT=CNAT +CNAP( I ) * ( DB ( I )/DB(NSEG) )**Z. 

200  CONTINUE 

WRITE(6,205)CNAT 

205  FORMAT ( 1 H  ,5X,'CNA  FOR  THE  BODY  = ' , F10 . 5) 

CALCULATE  DISTANCE  FROM  BASE  OF  BODY  TO  FRONT  FACE  OF  SEGMENT 

SUM3=0 . 0 
DO  215  1  =  1, NSEG 

DBFS(NSEG"-1-I)=SUM3+H(NSEG+1-I) 

SUM3=SUM3+H(  NSEG+l-I ) 

215  CONTINUE 

CAiCULATE  DISTANCE  FROM  SEGMENT  FRONT  FACE  TO  MOMENT  REFERENCE  AXIS 

IF ( RA . EQ  0 . 0)RA=CG0B 
DO  220  1  =  1. NSEG 
DBF  SR A ( I )=DBFS( I ) -R A 
220  CONTINUE 

CAl CUL A^E  CMA  FOR  EACH  SEGMENT( BASED  ON  AREA  AND  BASE  DIAMETER) 

DO  225  1  =  1  .NSEG 

CMA;  P=CMAP(  I  )«-DBFSRA(I)/DB(I)*CNAP(I) 

225  20NTINUE 

CALCULATE  CMA  f OR  EACH  SEGMENT ( BASED  ON  AREA  AND  SEGMENT  LENGTH) 

DO  51 C  I  =  i  , NSEG 

:mal  v : '=cmap • n*oe( i )/h( i )+dbfsra( i )/h( i )*cnap( i ) 

5 1 :  CONTINUE 

;a.:,,,a-e  CMA  tor  the  body(Based  on  area  and  base  diameter) 

:mat-o  : 

DC  2 32  1  =  1, NSEG 

CMAT  =CMAT -CMA( I ) * ( DBt I )/DB( NSEG) ) **3 . 


2  3"  lON’INJE 

A.  '  4  *t  THE 

LENGTH 

OP  THE  BODY 

0 

5v  . 

00  0  35  ; - 1 

NSEC. 

R ,  5 :  *  m  '  I 

2-£  '  ;n  f 

a  .a*;  HA 

-c  thE 

BODv( BASED  ON 

I 

■  ha  ■  ■:  o 

:  ■ : .  n  ■  e  g 


JNCl ASSIFIFP 


.v.v. 


v.v. 


uou  u  uu 


UNCLASSIFIED 


- 1 


CMATL=CMATL+CMAL( I)*(DB( I )/DB(NSEG))**2.*(H( I )/BL ) 

500  CONTINUE 

CALCULATE  THE  CENTER  OF  PRESSURE 

XM=(BL-RA)/BL 
XCP=XM-CMATL/CNAT 
WRITE  (6,240) 

240  FORMAT ( 1H  , 6X, 1 ISEG' ,8X, 'DBFS' ,10X, ' DBFSRA ' ,12X, ' CMA ' , 1 1 X , 

1  ' CMAL ' ,/) 

DO  245  1=1 ,NSEG 

WRITE (6 ,250) ISEG( I ) ,DBFS( I ) ,DBFSRA( I ) ,CMA( I ) , CMAL( I ) 

250  F0RMAT(1H  , 5X , 15 ,4( 5X ,F10. 5) ) 

245  CONTINUE 

WRITE(6 ,255)CNAT 

255  FORMAT ( 1H  ,5X, 'TOTAL  CNA  FOR  THE  BODY  =',F10.5) 

WRITE(6 ,260)CMAT 

260  F0RMAT(1H  ,5X,'T0TAL  CMA  FOR  THE  BODY  (BASED  ON  AREA  AND 
1  'BASE  DIAMETER  =' .F10.5) 

WRITE(6 , 505)CMATL 

505  F0RMAT(1H  , 5X , ' TOTAL  CMA  FOR  THE  BODY  (BASED  ON  AREA  AND 
1  'BODY  LENGTH  ='  ,F10.5) 

WRITE(6 ,265)BL 

265  F0RMAT(1H  ,5X, 'LENGTH  OF  THE  BODY  =',F10.5) 

WRITE(6,270)XCP 

270  FORMAT( 1H  ,5X, 'CENTER  OF  PRESSURE  AS  A  FRACTION  OF  BOOY  ', 

1  'LENGTH  =' ,F10.5) 

CALCULATE  GYROSCOPIC  STABILITY  AND  SPIN  RATE 

PI=3. 14159265 

IF( ISG.EQ.O)GO  TO  275 

W=(SG*PI*RHO*TIB*ABS(CMAT)*DB(NSEG)**3./(2.*AIB**2.))**.5*VEL 

RPM=W*30/PI 

WRITE(6 ,280)W,RPM 

280  FORMAT ( 1H  , 5X , 1  SPIN  RATE  (RAD/SEC)  =',F15.5,5X, 

1  'SPIN  RATE  (RPM)  =' .F15.5) 

GO  TO  285 

275  SG=2.*AIB**2.*(W/VEL)**2./(PI*RH0*TIB*CMAT*DB(NSEG)**3.) 
WRITE(6,290)SG 

290  FORMAT( 1H  , 5X , ' GYROSCOPIC  STABILITY  FACTOR  =',F10.5) 
WRITE(6,310)RH0,VEL 

310  FORMAT ( 1H  , 5X , 'AIR  DENSITY  =' ,F10 . 5 , 5X ,' PROJECTILE  VELOCITY  =' , 
1  F10.5) 

CALCULATE  COP  FOR  THE  FIRST  SEGMENT 

SPHERICAL  FIRST  SEGMENT 

.35  [Ff ITYPE .EQ. l)GO  TO  600 
CALL  F 4 231 66 ( BLUNT , COP ( 1) ) 

CO  TO  605 

si' A;  ( TRUNCATED)  FIRST  SEGMENT 


UNCLASSIFIED 


(“><“>(-)  o  o  o 


*»'  t. 


4 


.  A  -4 

.A, 


A  » 


'  A  A 

u  -s  ku  s ;  *»'.  * 


A  4  •  *  '•••('  A  U 

S'  Si  t 

>11  *  ■  *~ 

■  UHf,  '  t 


A  1  '  .  I  - .  •  S  S  )♦*»••  •• 

A  N  --  M  '  M 

A  >  4  .  1  .  AE  '  !•( 

A  A  '•  -!  •  V.  '  S  iMB( 


-i  -Mi  ,( 

Mt  (if 

;(.  :.Mf  .  ’  h( 

.Al_  1  4.61c*-  »E  .C* 
ut.  f  4  2  3 1  *  r  »E  '.E'TE 

C  CALCULATE  fRONTAL  ANC  rfEHEO  AREA' 

r 

?I  =  3  1416926 

S8=P:,0B(NlJt;,/*0BiNSEG;/4 

:  p r  i  type  rj  ;  >r,o  w  6 20 

AREA f  i  )=2  *P!*R2(  1  )*H(  1  j 
GO  TO  62 6 

62G  AREA(  1  )=PI*i  fRRc'fl  )-R2(  1  i  )**:  ♦Mfl)*H(l))**  6*  (  RR2U  )  *R2(  1  . 
626  DO  630  I=2,NSEG 

AREA( I ) =P I •( (RR2f ! )-R2( I ) J**2  5*( RR2( I J *R2(  I ,  , 

630  CONTINUE 

ss=o  0 

DO  636  1  =  1  ,NSEG 
SS=55+AREA(  I  ; 

636  CONTINUE 


CALCULATE  CDF,  BOOY  SKIN  FRICTION  DRAG  COEFFICIENT 


CDF=1 . 02*CF*CFCCF*SS/SB 

CALCULATE  CDO,  ZERO  LIFT  DRAG  COEFFICIENT 

SUMD=0.0 
DO  640  1=1 ,NSEG 

SUMD=SUMD+CDP( I ) * ( DB( I )/DB ( NSEG ) ) **2 . 

640  CONTINUE 

CALL  F423160( VEL ,CDB) 

CDO=CDF+SUMD+CDB 
C 

C  CALCULATE  CNQP .PITCHING  DERIVATIVE  FOR  SEGMENT 


UNCLASSIFIED 


o  o  o  nor 


JNCLASSIF  HP 


C-]‘ 


-A  -  JN  RASE  ART  A'.  ANL  BASf  D I AMI  T  t  R c  Rt  F Ft RPE  D  TO  FORWARD 

;  tr  I  ■  i;  1  t  (jMf  N  ’ 


S  NQR '  i 
CNQP  ]  ,  , 


a  A’f  'MQF  l  '  H 1 NG  Of  R  I  v  A '  !  v  (  FOR  SEGMENT 
-A  JL  'JN  BASf  AM  A'  ANL  SQUARE  Of  BASE  D I AMf  T (  R S 
RE  F  F  RRF  ;  *  rjRnrARL’  ‘  At.f  QE  THE  SEGMENT 

’  *  F  T  »  '  j  j’ ,  '  fj  tjtj  j 

a  i  3b  mqf  i 

jL  *t‘ 

MS  SA  F  n\,\  "A  'M  ’Af  ;  ,  ,AD(  1  ;  .CMQPM  ,  ; 

M;‘  Ds  R  / (j  NSfG 

'a  t  uin;  <  T«E  ’ai  i :  ad  <  i  i.cmqp.'  i  /  ;• 

*  7S  QNT ! NUf 

TRANSFER  thE  INDIVIDUAL  DERIVATIVES  TO  A  COMMON  MOMENT  CENTER 

DO  ATS  I  -1  .NSEC, 

R-08F$RA(  I  /DB( I  ; 

CMQ',  I  j=CMQP'  I  ,  -2  *R*CMAP(  I  j ♦R*CNQP(  I  )  -2  *R*R*CNAP(  I  ) 
f/S  rQNTINUE 

r  C A .  r , j L  a T f  CMQ  Ftjp  The  ENTIRE  800Y 

S  BASED  ON  A  COMMON  REFERENCE  AREA  AND  SQUARED  DIAMETER 

r 

CMQ T -0  C 
DO  680  1-1  . NSEC, 

CMQT-CMQT *CMQ(  }  )*(DB(  I  )/D8{ NSEC.) ) **4 
660  CONTINUE 

C  CALCULATE  THE  MAGNUS  MOMENT  COEFFICIENT 

C 

CALL  MAGNUS(BL .DB(NSEG) ,RES , CGOB , CMPA) 

C 

CALCULATE  CMAD.  THE  ACCELERATION  DERIVATIVE  CMAD  IN  THE  HYPERSONIC 
SPEED  RANGE  IS  EQUAL  TO  ZERO  WHEN  DETERMINED  BY  NEWTONIAN  THEORY 

CMAD=0 . 0 

CALCULATE  SO,  DYNAMIC  STABILITY  FACTOR 

AKB=(  AIB/( TMASS*D8(NSEG)**2 . )  )** .  5 
TKB=(TIB/(TMASS*DB(NSEG)**2. ) )** . 5 

SD=(2 . *(CNAT-CD0)+2 . *AKB**(-2 . )*CMPA)/(CNAT-2. *CD0-TKB**(-2 . )* 

1  (CMQT+CMAD) ) 


.  ■  «f .  ;  ,  .  j'.  *!.  t4 1 

A.  A  b>  JN '  ■  NQR  I  , 

,  ",  '  1 

•A-  A  .  -A  'A  .  AD'  I 

-  *.(.  ■  :  nsf  , 

A.  -A  7Hf  'A  :  AD?  ! 

-  1  v  '  IN  ’  :  N  li 


UNCLASSIFIED 


o  o  o 


UNCLASSIFIED 


l- 


C  WRITE  DYNAMIC  STABILITY  DATA 
C 

WRITE(6,700) 

700  F0RMAT(1H  , 10X , 1  *****  OYNAMIC  STABILITY  DATA  *****',//) 
WRIT£(6,705)ESR,AR 

705  FORMAT ( IH  , 5X EQUIVALENT  SAND  ROUGHNESS  =',E12.6,5X, 

1  'ADMISSIBLE  ROUGHNESS  = ' ,E 12 . 6) 

WRITE(6,710)RES,U 

710  F0RMAT(1H  , 5X REYNOLDS  NUMBER  =’,E12.6,5X, 

1  'ABSOLUTE  VISCOSITY  =',E12.6) 

WRITE(6,810)CRE 

810  F0RMAT(1H  ,5X, 'CUTOFF  REYNOLOS  NUMBER  =',E12.6) 
WRITE(6,715)SB,SS 

715  FORMAT( 1H  ,5X,'BOOY  MAXIMUM  FRONTAL  AREA  =',F10.5,5X, 

1  'BODY  WETTED  AREA  =' .F10.5) 

WRITE(6,720)AKB,TKB 

720  FORMAT( 1H  ,5X, 'DIMENSIONLESS  AXIAL  RADIUS  OF  GYRATION  =',F10.5, 
1  5X, 'DIMENSIONLESS  TRANSVERSE  RADIUS  OF  GYRATION  =',F10.5,//) 

DO  725  1=1 ,NSEG 
WRITE(6, 800)1 

800  FORMAT( 1H  ,5X, 'PROPERTIES  FOR  SEGMENT  NUMBER  ',15) 

WRITE(6 ,730)CDP( I ) ,CNQP( I ) 

730  FORMAT( 1H  ,5X,'CDP  =' ,F10.5,5X, 'CNQP  =',F10.5) 

WRITE(6,735)CMQP( I) ,CMQ(I) 

735  F0RMAT(1H  ,5X,'CMQP  =' ,F10.5,5X, 'CMQ  =',F10.5,//) 

725  CONTINUE 

WRITE(6,740)CF 

740  FORMAT( 1H  ,5X ,' INCOMPRESSIBLE  FLAT  PLATE  SKIN  FRICTION 
1  'COEFFICIENT  (CF)  =',F10.5) 

WR I TE ( 6 , 745 )CFCCF 

745  F0RMAT(1H  ,5X , 'COMPRESSIBLE/INCOMPRESSIBLE  SKIN  FRICTION  1 , 

1  'COEFFICIENT  (CFCCF)  =',F10.5) 

WRITE(6,750)CDF 

750  FORMAT( 1H  ,5X,'B0DY  SKIN  FRICTION  DRAG  COEFFICIENT  (CDF)  =', 

1  F10.5) 

WRITE(6,755)C0B 

755  FORMAT( 1H  ,5X,'B00Y  BASE  DRAG  COEFFICIENT  (CDB)  =',F10.5) 
WRITE(6, 760)C00 

760  F0RMAT(1H  ,5X,'B0DY  ZERO  LIFT  DRAG  COEFFICIENT  (CDO)  =',F10.5) 
WRITE(6,765)CMQT 

765  FORMAT( 1H  ,5X,'B0DY  PITCHING  DERIVATIVE  (CMQT)  =',F10.5) 
WRITE(6,780)CMPA 

780  F0RMAT(1H  ,5X, 'DERIVATIVE  OF  THE  MAGNUS  MOMENT  COEFFICIENT 
1  ' (CMPA)=  1 ,F10.5) 

WRITE(6,785)CMAD 

785  F0RMAT(1H  ,5X,'B0DY  ACCELERATION  DERIVATIVE  (CMPA)  =',F10.5) 
WRITE(6,790)SD 

790  F0RMAT(1H  ,5X,  'DYNAMIC  STABILITY  FACTOR  (SD)  =\F10.5) 

CALCULATE  STABILITY  OF  THE  BODY 

CALL  F611  (CMAT,SD,RH0,TIB,AIB,D8(NSEG),VEL,DEN) 

STOP 
END 


I 


UNCLASSIFIED 


O  C~)  o 
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SUBROUTINE  IHSS( ISEG.R2.R1 ,  H, DEN, V, AMASS, YB.B.AI ,TI) 

C 

C  CALCULATE  THE  INERTIA  PROPERTIES  FOR  A  HOLLOW  SPHERICAL  SEGMENT 
C 

C  ISEG  =  SEGMENT  NUMBER 
C  R2  =  OUTER  RAOIUS 
C  R1  =  INNER  RAOIUS 
C  H  =  HEIGHT  OF  SEGMENT 
C  DEN  =  DENSITY  OF  MATERIAL  (MASS/VOLUME) 

C  V  =  VOLUME 
C  AMASS  =  MASS 

C  YB  =  DISTANCE  FROM  REFERENCE  AXIS  TO  CENTROID 
C  B  =  DISTANCE  FROM  BASE  OF  SEGMENT  TO  CENTROID 
C  AI  =  AXIAL  (POLAR)  MOMENT  OF  INERTIA 
C  TI  =  TRANSVERSE  MOMENT  OF  INERTIA 
C 

PI=3 . 1415926 

V=PI*(R2*Rl*Rl+H'*R2*R2-H*Rl*Rl-R2**3./3.-2. *Rl**3./3. ) 

YB=3. *((R2*R2-R1'*R1)*(4. *H*R2-2. *H*H)-(R2*R2-R1*R1)**2. ) 

1  /(4 . *( (R2*R2-R1*R1 )*(3 . *H)+3 . *R2*R1*R1-R2**3 . -2 . *R1**3 . ) ) 
AMASS=DEN*V 
B=YB-R2+H 

AI=DEN*PI/2 .*(8./15.*(R2**5.-Rl**5. )+(R2-H)*(Rl**4 .-R2**4 . ) 

1  +2./3.*(R2-H)**3.*(R2*R2-Rl*Rl)) 

RR1=R2-R1 

RR2=(R2*R2-R1*R1) 

RR3=(R2**3.-R1**3. ) 

RR4=(R2*M.-R1**4.) 

RR5=(R2**5  — R 1**5  ) 

TI=DEN*PI*( ,25*RR4*(R1-R2+H)+RR2*(1./6.*(R1**3.-(R2-H)**3.) 

1  +YB*YB*(R1-R2+H)-YB*(R1*R1-(R2-H)**2 . )  )  +  .25*R2**4  .’“'RRl 

2  +R2*R2*(l./6.*RR3+YB*YB*RRl-YB*RR2)-3./20.*RR5+.5*YB*RR4 

3  -YB*YB/3 . *RR3) 

WRITE(6,300)ISEG 

300  F0RMAT(1H  ,5X,’  PROPERTIES  FOR  SEGMENT  NUMBER  ’,15) 
WRITE(6,301) 

301  F0RMAT(1H  , 10X , 1  HOLLOW  SPHERICAL  SEGMENT1) 

WRITE(6,302)R2,R1  ,H 

302  F0RMAT(1H  ,5X,'OUTER  RAOIUS  =  ' ,F10.5,5X, 1  INNER  RADIUS  = 

1  F10. 5, 5X, 'HEIGHT  =  \F10.5) 

WRITE(6,303)0EN,V, AMASS 

303  FORMAT( 1H  , 5X , ' DENSITY  =  ' ,F15. 8, 5X, 'VOLUME  =  ' ,F15.8,5X, 

1  'MASS  =  ' , F 1 5 .8) 

WRITE(6,304)YB 

304  FORMAT( 1H  , 5X, 'DISTANCE  FROM  REFERENCE  AXIS  TO  CENTROID  = 

1  F15.8) 

WRITE(6,305)B 

305  FORMAT! 1H  , 5X ,' DISTANCE  FROM  BASE  OF  SEGMENT  TO  CENTROID  = 

1  F15.8) 

WRITE(6,306)AI 

306  FORMAT! 1H  ,5X, 'AXIAL  MOMENT  OF  INERTIA  =  ' ,F15.8) 


UNCLASSIFIED 


o  r>  r> 


UNCLASSIFIED 


C  -  1ft 


WRITE(6,307)TI 

307  F0RMAT(1H  , 5X TRANSVERSE  MOMENT  OF  INERTIA  =  ' ,F1S  8) 
RETURN 
END 


SUBROUTINE  IHTC( ISEG.RR2.RR1 ,R2,R1 .H.OEN.V.AMASS.YB.AI ,TI ) 
C 

C  CALCULATE  THE  INERTIA  PROPERTIES  FOR  A  HOLLOW  TRUNCATED  CONE 
C 

C  ISEG  =  SEGMENT  NUMBER 
C  RR2  =  OUTER  RAOIUS  AT  LARGE  END 
C  RR1  =  INNER  RAOIUS  AT  LARGE  END 
C  R2  =  OUTER  RAOIUS  AT  SMALL  END 
C  R1  =  INNER  RADIUS  AT  SMALL  END 
C  H  =  HEIGHT  OF  TRUNCATED  CONE 
C  DEN  =  DENSITY  OF  MATERIAL  (MASS/VOLUME) 

C  V  =  VOLUME 
C  AMASS  =  MASS 

C  YB  =  DISTANCE  FROM  CONE  BASE  TO  CENTROID 
C  AI  =  AXIAL(POLAR)  MOMENT  OF  INERTIA 
C  TI  =  TRANSVERSE  MOMENT  OF  INERTIA 
C 

PI=3. 1415926 

DUMB=(RR2*RR2-RR1*RR1+R2*R2-R1*R1+RR2*R2-RR1*R1 ) 

V=PI*H/3. *OUMB 

YB= . 25*H*(RR2*RR2-RRl*RRl+3 . *R2*R2-3 . *Rl*Rl+2 . *RR2*R2- 
1  2. *RR1*R1)/DUMB 
AMASS=DEN*V 

AI=0EN*PI*H/10.*((RR2*M.+RR2**3.*R2+RR2*RR2*R2*R2+ 

1  RR2*R2**3.+R2**4.)-(RR1**4.+RR1**3.*R1+RR1*RR1*R1*R1+ 

2  RR1*R1**3.+R1**4. )) 

A1=RR1-R1 

A2=(RR1-R1)**2. 

A3=(RR1-R1)**3. 

A4=(RR1-R1)**4. 

B1=RR2-R2 

B2=(RR2-R2)**2. 

B3=(RR2-R2)**3 . 

B4=(RR2-R2)**4 . 

TI=( . 75*AMASS*(RR2**4 .-2 . *RR2**3 . *Bl+2 . *RR2**2 . *B2-RR2*B3 

1  +B4/5 . -RR1**4 . +2 . *RRI**3 . *Al-2 . *RR1**2 . *A2+RRl*A3-A4/5 . ) 

2  +3 .  *AMASS*(H*H*(RR2*RR2/3 .  +B2/5 .  -RR2*Bl/2 .  -RRl*RRl/3 . 

3  -A2/5 ,+RRl*Al/2 . )+Y8*YB*(RR2*RR2+B2/3.-RR2*Bl-RRl*RRl 

4  -A2/3 . +RRl*Al)+YB*H*(-RR2*RR2-B2/2 . +4 . *RR2*Bl/3 . 

5  +RRl*RRl+A2/2 .-4 . *RRl*Al/3 • ) ) )/DUMB 
WRITE(6 ,310) ISEG 

310  F0RMAT(1H  ,5X,'  PROPERTIES  FOR  SEGMENT  NUMBER  ',15) 

WRITE (6,311) 

311  F0RMAT(1H  , 10X , ' HOLLOW  TRUNCATED  CONE') 
WRITE(6,312)RR2,RR1 

312  F0RMAT(1H  ,5X,' OUTER  RADIUS  AT  LARGE  END  =  ',F10.5,5X, 

1  'INNER  RADIUS  AT  LARGE  END  =  \F10.5) 


UNCLASSIFIED 


ooooooo  r->  c->  o  o 
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W«  I  T  (  ;  6  .  U.|  )«  /  .HI 

ill  f  f)RMA  T  (  1 H  ,Sx.  UUTFR  MAD  I  US  At  SMAll  f  ND  ‘  .  F  10  S.SX, 

1  INNER  RADIUS  AT  SMAll  F ND  ‘  . F 1 0  Si 
WR I T  E  (  6  ,  314  )  m 

) i 4  FORMAT!  1M  .Sx.  HEIGHT  -  *  . F 1 0  S) 

WRI’E  i  f«.31S  )DEN.\.  .AMASS 

1 1 F0RMAT(1H  , Sx . ‘ DF  NSI T  Y  .FIS  8. SX, 'VOLUME  =  .FIS  8.SX. 

1  MASS  -  .FIS  H, 

WRITE  (6.31b )YB 

316  FORMAT; 1H  .SX,  DISTANCE  FROM  CONE  BASE  TO  CENTROID  =  '.FIS  8) 
WRITE(6.31/)AI 

II  ?  FORMAT! IN  .SX. 'AXIAL  MOMENT  OF  INERTIA  =  .FIS  8) 
WRITE(6.318)TI 

318  FORMAT; 1H  . Sx . ' TRANSVERSE  MOMENT  OF  INERTIA  =  '.FIS. 8) 

RE  TURN 

END 

C 

C . . . . . . 

c 

C***** . . 

c 

SUBROUTINE  C0EFF1 

GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.2.1.1-23 
NORMAL-FORCE-CURVE  SLOPE  FOR  SPHERICAL  SEGMENTS 

C0HM0N/XY1/X1 ( ll).Yl(ll) 

C0MM0N/BCD1/81 ( 11) ,C1( 11) ,D1( 11) 

DATA( Xl( I ). 1=1, 11 )/0. 0000, 0.1000, 0.2000 .0.3000, 0.4000, 

1  0 . 5000 , 0 . 6000 , 0 . 7000 , 0 . 8000 , 0 . 9000 , 

2  1 . 0000/ 

DATA( Yl( I ), 1=1, 11 )/0. 0000, 0.1748, 0.3497, 0.5070, 0.6399, 

1  0. 7517,0. 8392 ,0 . 9091 ,0.9580,0.9895, 

2  1 . 0000/ 

N=ll 

CALL  SPL I NE (N,X1 ,Y1 ,B1 ,C1 ,01) 

RETURN 
END 


SUBROUTINE  COEFF2 

GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.2.1.1-26 
NORMAL-FORCE-CURVE  SLOPE  FOR  CONE  FRUSTRUMS 

C0MM0N/XY2/X2( 11 ,10) ,Y2( 11 ,10) 

COMMON/BCD2/B2( 11, 10), C2 (11,10), D2 (11, 10) 

DIMENSION  B(10),C(10),D(10),XXX(10) ,YYY( 10) 

0ATA( X2( 1 , I ) , 1=1 , 10)/0. 0000 , 10 . 000 , 20 . 000 , 30 . 000,40 . 000 , 
1  50.000,60.000,70. 000 , 80 . 000 , 90 . 000/ 

0ATA( Y2( 1,1) ,1=1 , 10)/2 . 0000 , 1 . 9406 , 1 . 7727 , 1 . 5140 , 1 . 1818 , 
1  0.8287,0.5035,0.2448,0.0664,0.0000/ 

DATA( X2(2 , I ) , I -1 , 10 )/0 . 0000 , 10 . 000,20.000,30.000,40.000 , 
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1  SO  OOO . 60 . 000 . 70 . 000 , 80 . 000 , 90 . 000/ 

DATA( Y2( 2.1  )  .1  =  1,10)/!. 9825.1 .9161, 1.7448, 1.4790, 1.1678, 

1  0.8287.0.5035,0.2448.0. 0664 , 0 . 0000/ 

OATA( X2( 3. 1 ) ,  1  =  1 . 10)/0. 0000 . 10 . 000 . 20 . 000 , 30 . 000,40 . 000 . 

1  SO  000 , 60 . 000 , 70 . 000,80 . 000,90 . 000/ 

OATA( Y2( 3, I). 1=1, 10)/1. 9196, 1.8601, 1.6853, 1.4371, 1.1224, 

1  0.7867,0.4755,0.2448,0. 0664 , 0 . 0000/ 

OATA( X2( 4 , I ) , I =1 , 10 )/0 . 0000 . 10 . 000 , 20 . 000 , 30 . 000,40 . 000 , 

1  50 . 000 , 60 . 000 , 70 . 000 , 80 . 000 , 90 . 000/ 

OATA( Y2( 4 , I ) , I = 1 . 10 )/l . 821 7 , 1 . 7622 ,1.6119,1. 3601 , 1 . 0664 , 

1  0.7517,0.4510,0.2203,0.0629,0.0000/ 

DATA(X2( 5 , I ) , 1=1 , 10)/0. 0000 , 10.000,20 . 000,30.000,40.000, 

1  50 . 000 , 60 . 000 , 70 . 000 , 80 . 000 , 90 . 000/ 

DATA( Y2( 5, I), 1=1, 10)/1. 6853, 1.6329, 1.4755, 1.2587, 0.9825, 

1  0.6923,0.4231,0.1958,0.0594,0.0000/ 

DATA( X2( 6 , 1 ) , I  =  1 , 10)/0 . 0000 , 10 . 000 , 20 . 000 , 30 . 000,40 . 000 , 

1  50.000,60.000,70.000,80. 000 , 90 . 000/ 

0ATA( Y2(6, I). 1=1, 10)/1. 5070, 1.4510, 1.3217, 1.1259, 0.8811, 

1  0.6259,0.3706,0.1783,0.0559,0.0000/ 

DATA( X2( 7 , I ) , 1=1 , 10)/0 .0000 , 10 . 000,20.000,30 . 000,40.000 , 

1  50.000,60.000,70.000,80.000,90.000/ 

OATA( Y2( 7, I), 1=1, 10)/1. 2867, 1.2413, 1.1259, 0.9580, 0.7552, 

1  0.5315,0.3252,0.1538,0.0455,0.0000/ 

DATA( X2(8, i ) , 1=1 , 10)/0. 0000 , 10 . 000 , 20 . 000 , 30 . 000,40 . 000 , 

1  50.000,60.000,70.000,80.000,90.000/ 

DATA( Y2( 8 , I ) , 1=1 , 10 )/l . 0210,0. 9895 , 0 . 8986 , 0 . 7622 , 0 . 5909 , 

1  0.4161,0.2552,0.1224,0.0315,0.0000/ 

DATA( X2( 9 , I ) , 1=1 , 10 )/0 . 0000 , 10 . 000 , 20 . 000 , 30 . 000 , 40 . 000 , 

1  50.000,60.000,70. 000 , 80 . 000 , 90 . 000/ 

DATA(Y2(9, I ), 1=1, 10)/0. 7168, 0.6923, 0.6224, 0.5350, 0.4196, 

1  0.2937,0.1748,0.0804,0.0210,0.0000/ 

DATA(X2(10 , I) , 1=1 , 10 )/0. 0000, 10. 000 ,20 .000,30.000,40.000, 
1  50.000,60.000,70.000,80.000,90.000/ 

DATA(Y2( 10, I ), 1=1, 10)/0. 3671, 0.3566, 0.3287, 0.2797, 0.2203, 
1  0.1503,0.0909,0.0420,0.0105,0.0000/ 

DATA( X2 (11,1), 1=1, 10)/0. 0000 , 10 . 000 , 20 . 000,30 . 000,40 . 000 , 
1  50.000,60.000,70.000,80.000,90.000/ 

DATA(  Y2(  11 , 1 ) ,  1=1 , 10 )/0 . 0000 , 0 . 0000 , 0 . 0000 , 0 . 0000 , 0 . 0000 , 
1  0 . 0000 , 0 . 0000 , 0 . 0000 , 0 . 0000 , 0 . 0000/ 

N=10 


DO  200  J=1 ,11 
DO  210  K=1 , 10 
XXX(K)=X2(J,K) 

YYY(K)=Y2(J,K) 

210  CONTINUE 

CALL  SPLINE(N,XXX,YYY,B,C,D) 
DO  220  K=1 , 10 
B2(J,K)=B(K) 

C2(J ,K)=C(K) 

02(J ,K)=D(K) 

220  CONTINUE 
200  CONTINUE 
RETURN 
END 
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SUBROUTINE  COEFF3 

GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.2.2.1-25A 
CMAP  FOR  CONE  FRUSTRUMS 

COMMON/XY3/X3( 11 ,10),Y3(11,10) 

COMMON/XYB/XXB(6) ,YYB(6) 

COHHON/BCD3/B3( 11, 10), C3( 11, 10) ,03(11,10) 
COMMON/BCD3B/B3B(6) ,C3B(6) ,D3B(6) 

DIMENSION  B(10) ,C(10) ,0(10) ,XXX(10) , YYY( 10) 

DATA( X3(1,I),I=1,10)/15.279,20. 000,25 . 000,30 . 000,40 . 000, 

1  50.000,60.000,70.000,80.000,90.000/ 

0ATA( Y3( 1,1), 1=1, 10)/-2. 400, -1.800, -1.398, -1.149, -0.796, 

1  -0 . 564 , -0 . 388 , -0 . 249 ,-0.125,0. 0000/ 

DATA(X3(2, I ) , 1=1 , 10)/13 . 219,15. 000,20. 000,25. 000,30. 000, 

1  40.000,50.000,60.000,75.000,90.000/ 

DATA( Y3(2, I), 1=1, 10)/-2. 400, -2. 114, -1.599, -1.249, -1.017, 

1  -0 . 723 ,  -0 . 526 , -0 . 367,-0 .183,0. 0000/ 

DATA(X3( 3,1) ,I=1,10)/11 .159,15.000,20.000,25.000,30.000 , 

1  40.000,50.000,60.000,75.000,90.000/ 

DATA(Y3(3,I),I=l,10)/-2. 400, -1.807, -1.353, -1.083, -0.900, 

1  -0.658,-0.485,-0.350,-0.183,0.0000/ 

DATA(X3(4 , I ) , 1=1 , 10 )/9 .0990, 10.000,15 .000,20. 000,25.000, 

1  35.000,45.000,60.000,75.000,90.000/ 

DATA( Y3(4, I), 1=1, 10)/-2. 400, -2. 118, -1.453, -1.100, -0.903, 

1  -0.661,-0.516,-0.339,-0.170,0.0000/ 

DATA( X3( 5 , I ) , 1=1 , 10)/6 . 8670 , 10 . 000 , 15 . 000 , 20 . 000 , 25 . 000 , 

1  35.000,45.000,60.000,75.000,90.000/ 

DATA(Y3( 5, I), 1=1, 10)/-2. 400, -1.668, -1.159, -0.893, -0.737, 

1  -0.564,-0.460,-0.315,-0.159,0.0000/ 

DATA(X3(6, I ) , 1=1 , 10 )/5 .0000, 10.000 ,15.000,20.000,25.000, 

1  35.000,45.000,60.000,75.000,90.000/ 

DATA( Y3(6, I), 1=1, 10)/-2. 400, -1.239, -0.862, -0.685, -0.588, 

1  -0.471,-0.398,-0.287,-0.149,0. 0000/ 

0ATA(X3(7,I),I=l,10)/5. 0000, 10.000,15. 000, 20. 000, 30. 000, 

1  40.000,50.000,60.000,75.000,90.000/ 

DATA(Y3( 7, I), 1=1, 10)/-1. 595, -0.845, -0.613, -0.498, -0.402, 

1  -0.350,-0.308,-0.246,-0.138,0.0000/ 

DATA(X3(8, I ) , 1=1 , 10 )/5 .0000, 10.000, 15.000,20. 000 ,30.000, 

1  40.000,50.000,60.000,75.000,90.000/ 

DATA( Y3(8, I), 1=1, 10)/-0. 900, -0.519, -0.395, -0.339, -0.298, 

1  -0 . 284 ,-0.249,-0.204,-0.114,0. 0000/ 

DATA(X3(9 , I ) , 1=1 ,10)/5. 0000 , 10.000,15.000,25.000,35.000, 

1  45.000,55.000,65.000,75.000,90.000/ 

DATA( Y3( 9 , I ) , 1=1 , 10 )/-0 .436,-0.260,-0. 218,-0 .194,-0.190, 

1  -0.180,-0.163,-0.128,-0.083,0.0000/ 

DATA( X3( 10 , I ) , 1=1 , 10 )/5 .0000, 10.000,15 . 000,20 . 000,30.000, 
1  40.000,50.000,60.000,75.000,90.000/ 

DATA(Y3(10, I), 1=1, 10)/-0. 118, -0.090, -0.080, -0.083, -0.097, 
1  -0.100,-0.093,-0.080,-0.055,0.0000/ 

DATA(X3( 11 , I ) , 1=1 , 10)/0.0000 , 10.000,20.000,30.000,40.000, 
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1  50 . 000 , 60 . 000 , 70 . 000 , 80 . 000 , 90 . 000/ 

0ATA(Y3( 11, I), 1=1, 10)/0. 0000, 0.0000, 0.0000, 0.0000, 0.0000, 

1  0 . 0000 , 0 . 0000 , 0 . 0000 , 0 . OOOO , 0 . 0000/ 

DATA(XXB( I ) , 1=1 ,6)/5 . 0000 , 6 . 8670 , 9 . 0990 , 11 . 159 , 13 . 219 , 15 . 279/ 
DATA( YYB( I ) , 1=1 , 6 )/0 . 5000 , 0 . 4000 , 0 . 3000 , 0 . 2000 ,0.1000,0. 0000/ 
N=10 

DO  200  J=1 ,11 
DO  210  K=1 ,10 
XXX(K)=X3(J,K) 

YYY(K)=Y3(J,K) 

210  CONTINUE 

CALL  SPLINE(N,XXX,YYY,B,C,D) 

DO  220  K=1 , 10 
B3(J,K)=B(K) 

C3(J ,K)=C(K) 

D3(J,K)=D(K) 

220  CONTINUE 
200  CONTINUE 

GENERATE  THE  CUBIC  SPLINE  COEFFICIENTS  FOR  THE  BOTTOM  OF  THE 
GRAPH  WHERE  CMAP  =  -2.4  AND  A/D  IS  .5  TO  0 


CALL  SPLINE (N,XX8,YYB,B3B,C3B,D3B) 

RETURN 

END 

********************************************************** 

SUBROUTINE  COEFF4 

GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.2.2.1-25B 
CMAP  FOR  SPHERICAL  SEGMENTS 

COMMON/X Y4/X4 ( 1 3 ) , Y4 ( 1 3 ) 

COMMON/BCD4/B4 (13) ,C4(13) ,04(13) 

DATA(X4( I), 1=1, 13)/0. 0000, 0.0250, 0.0500, 0.1000, 0.2000, 

1  0.3000,0.4000,0.5000,0.6000,0.7000, 

2  0.8000,0.9000,1.0000/ 

DATA( Y4( I), 1=1, 13)/0. 0000, -.0918, -.1482, -.2188, -.3000, 

1  -.3600, -.3988, -.4341, -.4588, -.4800, 

2  -.4941, -.5012, -.5012/ 

N=13 

CALL  SPLINE(N ,X4 ,Y4 ,B4 ,C4 ,04) 

RETURN 


SUBROUTINE  F421123(BLUNT,CNAP) 


CALCULATES  CNAP  FOR  A  GIVEN  NOSE  BLUNTNESS 
NORMAL-FORCE-CURVE  SLOPE  FOR  SPHERICAL  SEGMENTS 


UNCLASSIFIED 
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C0MM0N/XY1/X1(11),Y1(11) 

C0MM0N/BCD1/B1 ( 11) ,C1(11) ,01(11) 

N=ll 

CALL  SEVAL( BLUNT, CNAP ,D1Y0,D2Y0,D3Y0,X1 ,Y1 ,B1 ,C1 ,D1 ,N) 

RETURN 

END 


SUBROUTINE  F421 126( THETA , CNAP ,AD) 

CALCULATES  CNAP  FOR  A  GIVEN  SEMICONE  ANGLE 
NORMAL-FORCE-CURVE  SLOPE  FOR  CONE  FRUSTRUMS 

C0MM0N/XY2/X2( 11 , 10),Y2(11,10) 

C0MM0N/BCD2/B2( 11 , 10),C2(11,10),D2(11,10) 

DIMENSION  CNA(ll) ,XXX(10) ,YYY(10) ,B(10),C(10),D(10) 
DIMENSION  XX(ll) ,BB(11) , CC( 11 ) , DD( 11 ) 

DATA(XX( I) , 1=1, ll)/0. 0,0. 1,0. 2,0. 3, 0.4, 0.5, 0.6, 0.7, 0.8, 
1  0.9, 1.0/ 

DO  200  J=1 ,11 
DO  210  K=1 , 10 
XXX(K)=X2(J,K) 

YYY(K)=Y2(J,K) 

210  CONTINUE 

DO  220  K=1 , 10 
B(K)=B2(J,K) 

C(K)=C2(J ,K) 

D(K)=D2(J,K) 

220  CONTINUE 

CALL  SE VAL( THETA, DUMB, 01Y0,D2Y0,D3Y0, XXX, YYY,B,C,D, 10) 
CNA(J)=DUMB 
200  CONTINUE 

CALL  SPLINE(11,XX,CNA,BB,CC,DD) 

CALL  SEVAL( AD , CNAP , D1 YO , D2Y0 , D3Y0 , XX , CNA , BB , CC , DD , 11 ) 

RETURN 

END 


A********************************************************* 


SUBROUTINE  F422125A(THETA,CMAP ,AD) 

CALCULATES  CMAP  FOR  A  GIVEN  SEMICONE  ANGLE 
CMAP  FOR  CONE  FRUSTRUMS 

COMMON/XY3/X3( 11 , 10) , Y3( 11,10) 

COMMON/XYB/XXB( 6 ) ,YYB(6) 

COMMON/BCD3/B3( 11 , 10 ) , C3 ( 1 1 , 10 ) , D3 ( 1 1 , 10 ) 
COMMON/BCD3B/B3B(6) ,C3B(6) ,D3B(6) 

DIMENSION  XXI (11) ,XX2(11) ,XX3(10),XX4(9),XX5(8),XX6(7) 
DIMENSION  XXX( 10) ,YYY( 10) ,B( 10) ,C( 10) ,D( 10) 

DIMENSION  CMA1( 11 ) ,BB 1(11), CC 1(11), DDl(ll) 

DIMENSION  CMA2(11) , BB2 (11) ,CC2(11) ,002(11) 

DIMENSION  CMA3( 10) ,BB3( 10) ,CC3( 10) ,DD3( 10) 
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DIMENSION  CMA4(9) ,BB4(9) ,CC4(9) ,DD4(9) 

DIMENSION  CMA5(8) ,BB5(8) ,CC5(8) ,DD5(8) 

DIMENSION  CMA6( 7) ,BB6(7) ,CC6( 7) ,DD6( 7) 

DATA(XX1( I) , 1=1 ,ll)/0.0,0. 1 ,0.2 ,0 . 3,0 .4 ,0 . 5,0 .6,0 . 7,0.8, 

1  0.9, 1.0/ 

REGION  1  (  THETA  .GE.  15.279) 

IF ( THETA. LT . 15.279)G0  TO  10 
DO  200  J=l,ll 
DO  210  K=1 , 10 
XXX(K)=X3( J ,K) 

YYY(K)=Y3( J ,K) 

210  CONTINUE 

DO  220  K=1 ,10 
B(K)=B3(J ,K) 

C(K)=C3(J,K) 

D(K)=D3( J ,K) 

220  CONTINUE 

CALL  SEVAL (THETA, DUMB ,01Y0 ,D2Y0 ,D3Y0 , XXX ,YYY,B ,C ,D , 10) 
CMA1(J)=DUMB 

200  CONTINUE 

CALL  SPLINE( 11 ,XX1 ,CMA1 ,BB1 ,CC1 ,DD1) 

CALL  SEVAL( AD , CMAP , D1Y0 ,02Y0 ,D3Y0 , XX 1 , CMA1 , BB1 . CC1 , 001 , 1 1 ) 
GO  TO  1000 

REGION  2  (  THETA  .LT.  15.279  BUT  .GE.  13.219) 

10  IF(THETA.LT. 13. 219)60  TO  20 
DO  201  J=2 ,11 
DO  211  K=1 , 10 
XXX(K)=X3( J ,K) 

YYY(K)=Y3( J ,K) 

211  CONTINUE 

DO  221  K=1 , 10 
B(K)=B3(J ,K) 

C(K)=C3(J,K) 

D(K)=D3(J,K) 

221  CONTINUE 

CALL  SEVAL ( THETA , DUMB , D 1 YO , D2Y0 , D3 YO , XX X , YYY , B , C , D , 1 0 ) 

CMA2(J)=DUMB 

XX2( J)=XX1( J) 

201  CONTINUE 

CALL  SEVAL( THETA .DUMB ,D1Y0 ,D2Y0 ,D3Y0 ,XXB ,YYB ,B3B ,C3B ,D3B ,6) 
CMA2 ( 1 ) =— 2 . 4 
XX2( 1 )=DUMB 

IF(AD.LT.DUMB)GO  TO  2000 

CALL  SPLINE ( 11 ,XX2 ,CMA2 ,BB2 ,CC2 ,DD2) 

CALL  SEVAL( AD ,CMAP ,01Y0 ,D2YO,D3YO , XX2 ,CMA2 ,BB2 ,CC2 ,DD2 , 11 ) 
GO  TO  1000 

REGION  3  (  THETA  .LT.  13.219  BUT  .GE.  11.159) 

20  IF (THETA.LT. 11 . 159 )G0  TO  30 
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DC  :.'4  MI.  10 
e ;  A  '  -  9  3  v  J  ,  A  ' 
k  -  - "  3  c  j  .  k  ' 

DiK  ■--33(J.O 
::4  cont:nue 

CA_  SEvAlv fHETA. DUMB, D1Y0,D2Y0,D3Y0, XXX, YYY,B,C,D,10) 

CMA5:  J- 3 ^ -DUMB 
»X5-  C  3^X*i(Jl 

204  CONTINUE 

CAl.  SEVAL (THETA, DUMB, D1Y0,D2Y0,D3Y0,XXB,YYB,B3B,C3B,D3B, 6) 

CMASi,  1  )=-2  4 
X X 5 .  I )=DUMB 

IF (AD. IT. DUMB )GO  TO  2000 

CALL  SPLINE(8,XX5,CMA5,BBS,CC5,DD5) 

CALL  SEVAL( AO ,CMAP ,D1Y0 ,D2Y0 ,D3Y0 ,XX5 ,CMA5 ,BB5 ,CC5 ,DD5 ,8) 

GO  TO  1000 

C  REGION  6  (THETA  . IT .  6.867  BUT  .GE.  5.000) 

SO  IF(IHETA.LT.5.000)GO  TO  60 
DO  205  J=o  ,  11 
DO  215  K- 1,10 
XXX(K)=X3(J,K) 

YYY(K)=Y3( J ,K) 

215  CONTINUE 

DO  225  K=1 , 10 
B(K)=B3(J,K) 

C(K)=C3(J ,K) 

D(K)=D3(J ,K) 

225  CONTINUE 

CALL  SE VAL( THETA, DUMB, D1 Y0,D2Y0,D3Y0, XXX, YYY,B,C,D, 10) 

CMA6( J-4)=DUMB 
XX6( J-4)=XX1( J) 

205  CONTINUE 

CALL  SEVAL( THETA, DUMB, DlY0,n2Y0,D3Y0,XXB,YYB,B3B,C3B,D3B, 6) 

CMA6(l)=-2.4 

XX6( 1)=DUMB 

IF(AD.LT.DUMB)GO  TO  2000 

CALL  SPLINE(7,XX6,CMA6,BB6,CC6,DD6) 

CALL  SEVAL(AD ,CMAP ,D1Y0 ,D2Y0 ,D3Y0 ,XX6 ,CMA6 ,BB6 ,CC6 ,DD6 ,7) 

GO  TO  1000 
60  CMAP=0 . 0 

WRITE (6, 4000) THETA 

4000  F0RMAT(1H  ,5X,'THETA  =  1  ,F10.5,5X, 'OUT  OF  THETA  RANGE  '  'J 
GO  TO  1000 

2000  WRITE (6, 3000) AD, THETA 

3000  FORMAT( 1H  ,5X,'A/0  =  ' ,F10 . 5 , 5X , 1  WHICH  IS  OUT  OF  PANG; 

1  '  FOR  THETA  =  ' ,F10.5) 

CMAP=999 . 99 
1000  RETURN 
END 
C 
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SUBROUTINE  F422125B(BLUNT,CMAP) 

CALCULATES  CMAP  FOR  A  GIVEN  NOSE  BLUNTNESS 
CMAP  FOR  SPHERICAL  SEGMENTS 

C0MM0N/XY4/X4(13) ,Y4(13) 

C0MM0N/BCD4/B4(13) ,C4(13) ,D4(13) 

N=13 

CALL  SE VAL (BLUNT , CMAP , 01Y0 , D2 YO , D3Y0 ,X4,Y4,B4,C4,D4,N) 

RETURN 

END 


********************************************************** 


SUBROUTINE  F423166(BLUNT,CDP) 

CALCULATE  CDP,  NEWTONIAN  DRAG  COEFFICIENT  FOR  SPHERICAL 
SEGMENTS  REFERRED  TO  BASE  AREA  OF  SEGMENT 

CDP=BLUNT**2 . -2 . *BLUNT+2 . 

RETURN 

END 

********************************************************** 


SUBROUTINE  F721213B(CMQP) 

CALCULATE  CMQP,  PITCHING  DERIVATIVE  FOR  SPHERICAL  SEGMENTS 
BASED  ON  THE  BASE  AREA  AND  THE  SQUARE  OF  THE  BASE  DIAMETER 


CMQP=-.5 

RETURN 

END 


********************************************************** 


SUBROUTINE  F721212(THETA,AD,CMQP) 

CALCULATE  THE  PITCHING  DERIVATIVE  CMQP  DUE  TO  INCLINED  SIDES 
OF  CONE  FRUSTRUMS 

PI=3. 1415926 
RAD=THETA*P 1/180. 

A=6 . *AD*AD*( 1 . -AD*AD) 

B=-8 . *AD*( 1 . -AD**3 . ) 

C=3 .  *(  1 . -AD**4 . ) 

CMQP=-1 ./(6 .*(SIN(RAD) )**2 . )*( A*(COS(RAD))**4.+B*(COS(RAD))**2. 
1  +C) 

IF(THETA.GE.5.0)GO  TO  100 

CMQP=0.0 

WRITE(6,105) 

105  F0RMAT( 1H  ,5X,'0UT  OF  THETA  RANGE') 


UNCLASSIFIED 


o  o  o  o  ooo  ooo  ooo  o  o  o  o  ooo  o  o  o 


UNCLASSIFIED 
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100  IF(CMQP.GE.-4.8)G0  TO  200 
WRITE(6,110) 

110  F0RMAT(1H  ,5X,'0UT  OF  CMQP  RANGE') 

200  RETURN 
END 
C 

Q***A**AA*4r*A;AA*AA4nlcA***AAA*A*ilr***AA***AA***AAAAr*4rA:**4nkA:**ir 

SUBROUTINE  F721213A(THETA , AD , CMQP ) 

CALCULATE  THE  PITCHING  DERIVATIVE  CMQP  FOR  THE  TOTAL  CONE  FRUSTRUM 

CALL  F721212( THETA, AD, CMQP) 

IF(CMQP.EQ.0.0)GO  TO  100 
CMQP=CMQP-(AD**4 . )/2 . 

100  RETURN 
END 


********************************************************** 


SUBROUTINE  F423167S(THETA,AD,CDP) 

CALCULATE  CDP,  DRAG-FORCE  COEFFICIENT  DUE  ONLY  TO  THE  INCLINED 
SIDES  OF  A  CONE  FRUSTRUM,  BASED  ON  BODY  BASE  AREA 

PI=3. 1415926 
RAD=THETA*P 1/180. 

CDP=2.*(SIN(RAD))**2.*(1 .-AD*AD) 

RETURN 

END 

A****************-*************************************** 


SUBROUTINE  F423167T(THETA,AD,CDP) 

CALCULATE  CDP,  DRAG-FORCE  COEFFICIENT  FOR  THE  TOTAL  CONE  FRUSTRUM 

PI=3. 1415926 
RAD=THETA*P 1/180. 

CDP=2. *((SIN(RAD))**2. *(1 .-AD*AD)+AD*AD) 

RETURN 

END 


********************************************************** 


SUBROUTINE  COEFF5 

GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4. 2. 3. 1.-68 
COMPRESSIBILITY  EFFECT  ON  TURBULENT  SKIN  FRICTION 

COMMON/XY5/X5(8),Y5(8) 

COMMON/BCD5/B5(8) ,C5(8) ,05(8) 

DATA( X5(I) , 1=1 ,8)/0. 0000 , 1 . 0000 , 2 . 0000 , 3 . 0000 , 

1  4.0000,5.0000,6.0000,10.000/ 


UNCLASSIFIED 


ooo  ooo  oooo  ooo  ooooo  oooo  ooo 


UNCLASSIFIED 


C-29 


DATA(Y5( I), 1=1, 8)/l. 0000, 0.9271, 0.7708, 0.6146, 
1  0.5000,0.3958,0.3333,0.1289/ 

CALL  SPLINE(8,X5,Y5,B5,C5,D5) 

RETURN 

END 


*********************************!![*********************** 


SUBROUTINE  F423168(VEL,CFCCF) 

CALCULATES  CFC/CF,  COMPRESSIBILITY  EFFECT  ON  TURBULENT 
SKIN  FRICTION 

C0MM0N/XY5/X5(8) ,Y5(8) 

C0MM0N/BCD5/B5(8) ,C5(8) ,05(8) 

VEL  =  PROJECTILE  VEL0CITY( IN/SEC) 

A  =  SPEED  OF  SOUND  IN  AIR(1116.4  FT/SEC) 

AMACH  =  PROJECTILE  VELOCITY  (MACH  NUMBER) 

A=1116.4 

AMACH=VEL/( 12 . *A) 

IF( AMACH. GT. 10. )G0  TO  100 

CALL  SEVAL ( AMACH , CFCCF , D1Y0 , D2Y0 , D3Y0 , X5 , Y5 , B5 , C5 , D5 , 8) 

GO  TO  200 

100  WRITE (6, 300) AMACH 

300  F0RMAT(1H  ,10X,'MACH  NUMBER  =' ,F10.5,5X, 'WHICH  IS  OUT  OF  RANGE') 
200  RETURN 
END 


***************  ****************************************** 


SUBROUTINE  C0EFF6 

GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.1.5.1-27 
CUTOFF  REYNOLDS  NUMBER 

C0MM0N/XY6/X6(5) ,Y6(5) 

COMMON/BCD6/B6( 5) ,C6(5) ,D6(5) 

DATA( X6( I ) , 1=1 , 5 )/l . 80E03 , 1 . 65E04 , 1 . 47E05 , 1 . 30E06 , 1 . 00E07/ 
DATA( Y6( I ) , 1=1 , 5)/l . 00E05 , 1 . 00E06 , 1 . 00E07 , 1 . 00E08 ,8 . 50E08/ 
CALL  SPLINE (5,X6 ,Y6,86,C6 ,D6) 

RETURN 

END 

********************************************************** 


SUBROUTINE  F415127(AR ,CRE) 

CALCULATES  CUTOFF  REYNOLDS  NUMBER 

C0MM0N/XY6/X6(5) ,Y6(5) 
COMMON/BCD6/B6(5) ,C6(5) ,D6(5) 


UNCLASSIFIED 


UNCLASSIFIED 


C-30 


C  AR  =  ADMISSIBLE  ROUGHNESS 
C  CRE  =  CUTOFF  REYNOLDS  NUMBER 
C 

CALL  SEVAL( AR , CRE , D1Y0 , 02Y0 , D3Y0 , X6 , Y6 , B6 , C6 , D6 , 5 ) 
RETURN 
END 
C 

Q********************************************************** 

C 

SUBROUTINE  COEFF7 


C 

C 

C 

C 

C 


GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.1.5.1-26 
TURBULENT  MEAN  SKIN  FRICTION  COEFFICIENT  ON  AN  INSULATED 
FLAT  PLATE 


1 

2 

3 

4 

5 

1 

2 

3 

4 

5 


COMMON/XY7/X7( 26 ) ,Y7(26) 

COMMON/BCD7/B7 ( 26 ) ,C7( 26) , D7( 26) 

DATA(X7(I) , I=1,26)/3.3E05,5.0E05,7.0E05,9.0E05,1.0E06, 

1 . 5E06 , 2 . OE06 , 3 . 0E06 , 5 . 0E06 , 7 . 0E06 , 
9 . 0E06 , 1 . 0E07 , 1 . 5E07 , 2 . 0E07 , 3 . 0E07 , 
5 . 0E07 , 7 . 0E07 , 9 . 0E07 , 1 . 0E08 , 1 . 5E08 , 
2 . 0E08 , 3 . OE08 , 5 . 0E08 , 7 . 0E08 , 9 . 0E08 , 
1.0E09/ 

DATA(Y7( I ) , 1=1 , 26)/ . 00550 , . 00505 , . 00474 , . 00454 , . 00445 , 

.00415, .00395, .00369, .00338, .00320, 
. 00306 , . 00300 , . 00283 , . 00270 , . 00254 , 
. 00235 , . 00223 ,. 00215 , . 00212 , . 00199 , 
. 00191 , . 00182 , . 00170 , . 00164 , . 00160 , 
.00158/ 

CALL  SPLINE(26,X7,Y7,B7,C7,D7) 

RETURN 

END 


Q*  *  irk  ★  *  ★  **  ft  ft  it  *  *  ★  *  ★  *  ★  ft  *  ★  *  ★  ★  **  *  *  *  *  *  *  ★  *  *  *  it  Hr  Hr  ft  ★  ★  ★  it  It  It  it  it  *  *  ★  *  irk 

C 

SUBROUTINE  F415126(RE ,CF) 

C 

C  CALCULATE  TURBULENT  MEAN  SKIN  FRICTION  COEFFICIENT  ON  AN 
C  INSULATED  FLAT  PLATE 
C 

COMMON/XY7/X7 ( 26 ) ,Y7(26) 

COMMON/BCD7/B7(26) , C7 ( 26 ) , D7( 26 ) 

C 

C  RE  =  REYNOLDS  NUMBER 

C  CF  =  SKIN  FRICTION  COEFFICIENT 

C 

CALL  SEVAL(RE,CF,D1Y0,D2Y0,D3Y0,X7,Y7,B7,C7,D7,26) 

RETURN 

END 


C 

C 

SUBROUTINE  C0EFF8 
C 


UNCLASSIFIED 


•mm 


ooo  ooo  o  o  o  o  o  o  ooo  ooo 


UNCLASSIFIED 
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C  GENERATES  CUBIC  SPLINE  COEFFICIENTS  FOR  FIG.  4.2.3.1-60 
C  BASE  DRAG  COEFFICIENT  WITH  NO  BOATTAIL 
C 

C0MM0N/XY8/X8(8),Y8(8) 

C0MM0N/BCD8/B8(8) ,C8(8) ,D8(8) 

0ATA(X8( I ) , 1=1 , 8)/l . 5000 , 2 . 0000 , 2 . 5000 , 3 . 5000 , 

1  4.5000,6.0000,7.0000,10.000/ 

DATA(Y8( I ) , 1=1 , 8)/0 . 1776 , 0 . 1440 ,0.1168,0. 0800 , 

1  0.0560,0.0352,0.0272,0.0110/ 

CALL  SPLINE(8 ,X8,Y8,B8,C8,D8) 

RETURN 

END 

********************************************************** 


SUBROUTINE  F423160( VEL ,CDB) 

CALCULATES  CDB,  BASE  DRAG  C0EFFICIENT(WITH  NO  BOATTAIL) 

COMMON  /XY8/X8(8) ,Y8(8) 

COMMON  /BCD8/B8(8) ,C8(8) ,08(8) 

VEL  =  PROJECTILE  VELOCITY  (IN/SEC) 

A  =  SPEED  OF  SOUND  IN  AIR  (1116.4  FT/SEC) 

AMACH  =  PROJECTILE  VELOCITY  (MACH  NUMBER) 

CDB  =  BASE  DRAG  COEFFICIENT 


A=1116.4 

AMACH=VEL/ ( 12 . *A) 

CALL  SEVAL ( AMACH , CDB , D1Y0 , D2Y0 , D3Y0 , X8 , Y8 , B8 , C8 , D8 , 8) 

RETURN 

END 


********************************************************** 


SUBROUTINE  F72119A( THETA, AD, CNQP) 

CALCULATE  CNQP,  PITCHING  DERIVATIVE  FOR  CONE  FRUSTRUMS 

PI=3. 1415926 
RAD=THETA*P 1/180. 

CNQP=2 ./(3.*TAN(RAD))*(2 .*(1-AD**3. )-3.*AD*(C0S(RAD)**2. ) 
1  *( 1-AD*AD) ) 

IF(THETA.GE.5.0)GO  TO  100 
CNQP=0.0 
WRITE(6,105) 

105  F0RMAT(1H  ,5X,'OUT  OF  THETA  RANGE') 

100  I F ( CNQP . LE .4 .4)G0  TO  200 
WR ITE (6,110) 

110  F0RMAT(1H  ,5X,'OUT  OF  CNQP  RANGE') 

200  RETURN 
END 
C 

Q******************************************* *************** 


UNCLASSIFIED 


ooo  oooo  ooo  ooo  ooo 


UNCLASSIFIED 


C 

SUBROUTINE  F72119B(BLUNT,CNQP) 

CALCULATE  CNQP,  PITCHING  DERIVATIVE  FOR  SPHERICAL  SEGMENTS 

CNQP=( BLUNT*( 2 . -BLUNT) ) ** . 5 
RETURN 
END 


***************************************** **********  ******* 


SUBROUTINE  MAGNUS(BL,DB,RES ,CG0B ,CMPA) 

CALCULATE  THE  DERIVATIVE  OF  THE  MAGNUS  MOMENT  COEFFICIENT 

DMCP=AB$( CGOB- . 333*BL ) 

CMPA=26 . 3* ( BL/DB ) **2 .  *0MCP/( RES** . 5 ) 

RETURN 

END 


********************************************************** 


SUBROUTINE  F611  (CMAT,SD,RHO,TIB,AIB,DB,VEL,DEN) 

SUBROUTINE  TO  DETERMINE  BODY  STABILITY 

DIMENSION  TLS( 25 ) ,HLS(25) ,RCS(25) ,AMASSL(25) 

PI=3. 1415926 

IF(SD.LT.O.O.OR.SO.GT.2.0)GO  TO  10 
IF(CMAT.GT.O.O)GO  TO  20 
WRITE(6,25)CMAT 

25  F0RMAT( 1H  ,5X,'CMA  =' ,F10. 5, 5X, ' (NEGATIVE) ' ) 

WRITE(6,30)SD 

30  F0RMAT(1H  ,5X,'SD  =' ,F10.5,5X, '(IN  THE  RANGE  0-2)') 
WRITE(6,35) 

35  F0RMAT( 1H  ,5X,'BODY  IS  DYNAMICALLY  STABLE  AT  ANY  SPIN  RATE') 
GO  TO  500 
20  SGI=S0*(2.-SD) 

SG=1 ./SGI 

W=(SG*PI*RHO*TIB*CMAT*DB**3 ./(2 . *AIB**2. ) )** . 5*VEL 
RPM=W*30./PI 
WRITE(6,40)CMAT 

40  F0RMAT(1H  ,5X,'CMA  =' ,F10. 5, 5X (POSITIVE) ' ) 

WRITE(6,30)SD 
WRITE(6,65)SGI ,SG 

65  F0RMAT( 1H  ,5X,'SGI  =' .F10.5.5X, 1 SG  =' ,F10. 5,5X , 1 (VALUES  ON  ' 
1  'DYNAMIC  STABLE/UNSTABLE  CURVE)') 

WRITE(6 ,45)W,RPM 

45  F0RMAT(1H  ,5X,'B0DY  IS  DYNAMICALLY  STABLE  AT  SPIN  RATES  ', 

1  'GREATER  THAN  ' .E12.6.3X, '(RAD/SEC) ' , 5X,E12.6,3X, ' (RPM) ' ) 

C 

C  CALCULATE  THE  EQUIVALENT  SPIN  RATE  OF  THE  SFF  LINER 
C  BASED  ON  THE  CONSERVATION  OF  ANGULAR  MOMENTUM 


UNCLASSIFIED 


o  o  o  <-> 


UNCLASSIFIED 


C-33 


READ( 5 , 100) ISP  IN 

100  F0RMAT( 15) 

IF( ISP  IN . NE .l)GO  TO  500 
WRITE(6,200) 

200  F0RMAT(1H  ,//, 15X LINER  DATA1,/) 

WRITE(6 ,205) 

205  F0RMAT(1H  ,8X, 1 TLS 1 ,13X, ' HLS 1 ,12X, 1 RCS 1 ,10X, 'AMASSL' ,5X, 

1  'SEGMENT',/) 

TAM=AIB*W 
READ(5,100)NLS 
TMASSL=0.0 
SUM1=0 .0 
DO  105  1=1, NLS 

READ(5,110)TLS(I) ,HLS(I) ,RCS(I) 

110  FORMAT(3F10.5) 

AMASSL( I )=DEN*TLS( I )*HLS( I)*RCS( I)*2. *PI 
TMASSL=TMASSL+AMASSL( I ) 

WRITE (6,210)TLS(I) ,HLS( I) ,RCS( I) ,AMASSL(I ) ,1 

210  F0RMAT(1H  ,3(5X ,F10 . 5) ,E12. 5, 15) 

SUM1=SUM1+AMASSL( I )*RCS( I )**2 . 

105  CONTINUE 
WL=TAM/SUM1 
RPML=WL*30./PI 
WRITE(6 ,215)TMASSL 

215  F0RMAT(1H  ,5X, 'TOTAL  MASS  OF  THE  LINER  =  ' ,E15.6,/) 
WRITE(6,220)WL ,RPML 

220  F0RMAT(1H  ,5X,'SFF  IS  DYNAMICALLY  STABLE  AT  CHARGE  SPIN  ', 

1  'RATES  GREATER  THAN  ' ,E12.6,3X, ' (RAD/SEC) ' ,E12.6,3X, '(RPM) 1 ) 
GO  TO  500 

10  IF(CMAT.LT.O.O)GO  TO  50 
WRITE(6,40)CMAT 
WRITE(6,55)SD 

55  F0RMAT(1H  ,5X,'SD  =' ,F10 . 5,5X, ' (NOT  IN  THE  RANGE  0-2)') 
WRITE(6 ,60) 

60  F0RMAT(1H  ,5X,'80DY  IS  DYNAMICALLY  UNSTABLE  AT  ANY  SPINRATE') 
GO  TO  500 

50  SGI=SD*(2.-SD) 

SG=1 ./SGI 
SGA=ABS(SG) 

W=(SGA*PI*RH0*TIB*ABS(CMAT)*DB**3 ./( 2 . *AIB**2 .))**. 5*VEL 

RPM=W*30./PI 

WRITE(6 ,25)CMAT 

WRITE(6,55)SD 

WR I TE (6,65) SG I ,SG 

WRITE(6,70)W,RPM 

70  F0RMAT(1H  ,5X,'B0DY  IS  DYNAMICALLY  STABLE  AT  SPIN  RATES  ', 

1  'LESS  THAN  ' , E 12 . 6 , 3X , ' (RAD/SEC) ' ,E12.6,3X, ' (RPM) ' ) 

500  RETURN 
END 


********************************************************** 

********************************************************** 


UNCLASSIFIED 


o  o  o 


r 


UNCLASSIFIED 


SUBROUTINE  SPLINE  (N,X,Y,B,C,D) 

C 

C  GENERATES  CUBIC  SPLINE  COEFFICIENTS 
C 

REAL  X(101)  ,  Y(101)  ,  B(101)  ,  C(101)  ,  D(101) 

NM1  =  N  -  1 
IF(N-2)60,50,1 

I  D(l)  =  X(2)  -  X(l) 

C( 2)  =  (Y(2)  -  Y(l))/D(l) 

DO  10  1=2, NM1 

0(1)  =  X( 1+1)  -  X(I) 

B(I)  =  2.*(D(I-1)  +  D( I ) ) 

C(I+1)  =  (Y(I+1)  -  Y(I))/D(I) 

C(I)  =  C(I+1)  -  C(I) 

10  CONTINUE 
B(l)  =  -0(1) 

B(N)  =  -D(N-l) 

C(l)  =  0.0 
C(N)  =  0.0 
IF(N-3)11,15,11 

II  C(l)  =  C(3)/(X(4)-X(2))  -  C(2)/( X( 3)— X( 1 ) ) 

C(N)  =  C(N-l)/(X(N)-X(N-2))  -  C(N-2)/(X(N-l)-X(N-3)) 

C(l)  =  C(1)*D(1)*D(1)/(X(4)-X(1)) 

C(N)  =  — C(  N )  *D(  N— 1 )  ’*rD(  N— 1 )/( X  ( N )— X(  N— 3 ) ) 

15  DO  20  1=2, N 

T  =  D(I-1)/B(I-1) 

B( I)  =  B(I)  -  T*D(I-1) 

C(I)  =  C(I)  -T*C(I-1) 

20  CONTINUE 

C(N)  =  C(N)/B(N) 

DO  30  IB  =  l.NMl 
I  =  N  -  IB 

C(I)  =  (C(I)  -  D(I)*C(I+1))/B(I) 

30  CONTINUE 

B(N)  =  (Y(N)  -  Y(NM1) )/D(NMl)  +  D(NM1)*(C(NM1)  +  2.*C(N)) 
DO  40  I  =  l.NMl 

B( I)  =  (Y( 1+1)  -  Y(I))/D(I)  -  D(I)*(C(I+1)  ♦  2.*C(I)) 
D(I)  =  (C( 1+1)  -  C(I))/D(I) 

C( I)  =  3.*C(I) 

40  CONTINUE 

C(N)  =  3.*C(N) 

D(N)  =  D(N-'i) 

RETURN 

50  B(l)  =  (Y(2)-Y(1))/(X(2)-X(1)) 

C(l)  =  0.0 
0(1)  =  0.0 
B(2)  =  B(l) 

C(2)  =  0.0 
0(2)  =  0.0 
60  RETURN 
END 


********************************************** ************ 


UNCLASSIFIED 


C-34 


o  o  o 


UNCLASSIFIED 


C-35 


SUBROUTINE  SEVAL(X0,Y0,D1Y0,D2Y0,D3Y0,  X,Y,B,C,D,N) 

EVALUATES  CUBIC  SPLINE  FOR  A  GIVEN  ABSCISSA 

REAL  X(100)  ,  Y(IOO)  ,  B(IOO)  ,  C(IOO)  ,  D(IOO) 

DATA  1/1/ 

IF( I-N)2 ,1 ,1 
11  =  1 

2  IF( XO— X( I ) )  10,3,3 

3  IF(X0-X(I+1))  30,30,10 
10  1  =  1 

J  =  N  +  1 

20  K  =  ( I+J)/2 

I F( XO— X( K ) )  21,22,22 

21  J  =  K 

GO  TO  23 

22  I  =  K 

23  I F( J  -  (I+l))30,30,20 
30  DX  =  XO  -  X( I) 

YO  =  Y(I)  +  DX*(B(I)  +  DX*(C(I)  +  DX*D(I))) 

D1Y0  =  B(I)  +  DX*(2.*C(I)  +  3.*D(I)*DX) 

D2Y0  =  2.*C(I)  +  6.*D(I)*DX 
D3Y0  =  6. *0(1) 

RETURN 

END 


UNCLASSIFIED 


UNCLASSIFIED 


C-36 


SFFSTAB  INPUT 


TITLE  CARD  ( 13A6,  2A1) 


BODY  MATERIAL  DENSITY  (F10.5) 


SEGMENT  CARD  (215) 


ITYPE  NSEG 


ITYPE  =  0  -  spherical  segments  and  truncated  cones  are  used 
1  -  only  truncated  cones  are  used 
NSEG  ■  number  of  segments 

INDIVIDUAL  SEGMENT  CARDS  -  AS  REQUIRED 

SPHERICAL  SEGMENT  CARD  (15,  5X,  3F10.5)  -  AS  REQUIRED 


ISEG  =  segment  number 
R2  =  outer  radius 
R1  *  inner  radius 
H  =  segment  height 

TRUNCATED  CONE  SEGMENT  CARD  (15,  5X,  5F10.5)  -  AS  REQUIRED 


RR2 

RR1 

R2 

R1 

H 

segment  number 
outer  radius  at  large  end 
inner  radius  at  large  end 
outer  radius  at  small  end 
inner  radius  at  small  end 
segment  height 

UNCLASSIFIED 


I 


UNCLASSIFIED 


C-37 


REFERENCE  AXIS  (F10.5)  RA  =  0.0  if  C  of  G  is  to  be  used 


10 


GYROSCOPIC  STABILITY  CARD  (15) 


ISG  =  0  -  calculate  SG 
=  1  -  calculate  W 


CARD  FOR  ISG  =  0  (3FI0.5)  -  AS  REQUIRED 


w 

11 

RHO 

D  2 

VEL 

0  3 

w//////////////////M/mm 

W  =  spin  rate  (RAD/SEC) 

RHO  =  air  density  (lb  sec2/in4) 
YEL  =  body  velocity  (in/sec) 


CARD  FOR  ISG  =  1  (3F10.5)  -  AS  REQUIRED 


SG 

RHO 

VEL 

M/M/MMIMMiim 

SG  =  gyroscopic  stability  factor 
DRAG  PARAMETERS  CARD  (2E10.5) 


UNCLASSIFIED 


C-38 


f 


(15) 


5 


I  SPIN  =  1  *  calculates  charge  spin  rate 
(15) 


NLS  =  number  of  liner  segments 
(3F10.5) 


ONLY 
—  FOR 

ISPIN  =  1 


Repeat  for  each  segment 

TLS  =  thickness  of  segment 
HLS  =  height  of  segment 
RCS  =  radius  of  centroid  of  crossectional 
area  of  liner  segment 


UNCLASSIFIED 


UNCLASSIFIED 


C-39 


SAMPLE  PROBLEM  INPUT 


BASELINE  DESIGN  SHOT  #01 
.000837 
0  5 


1 

.3564 

0.0 

.2152 

2 

.4048 

0.0 

.3193 

0.0 

.2024 

3 

.4576 

0.0 

.4048 

0.0 

.3107 

4 

.4704 

0.0 

.4576 

0.0 

1.1717 

5 

.5488 

0.0 

.4704 

0.0 

.6942 

0.0 


1 

1.0.000000115  92262. 

2.5932E-09  .16E-03 

1 

20 


.121 

.1 

.05 

.122 

.1 

.15 

.123 

.1 

.25 

.124 

.1 

.35 

.125 

.1 

.45 

.126 

.1 

.55 

.127 

.1 

.65 

.128 

.1 
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